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Recent experimental data and progress in nuclear structure modeling have lead to improved de- 
scriptions of astrophysically important weak-interaction processes. The review discusses these 
advances and their applications to hydrostatic solar and stellar burning, to the slow and rapid 
neutron-capture processes, to neutrino nucleosynthesis, and to explosive hydrogen burning. Spe- 
cial emphasis is given to the weak-interaction processes associated with core-collapse supernovac. 
Despite some significant progress, important improvements are still warranted. Such improve- 
ments are expected to come from future radioactive ion-beam facilities. 
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The weak interaction is one the four fundamental 
forces in nature. Like the other three - strong, elec- 
tromagnetic and gravitation - it plays a keyrole in many 
astrophysical processes. This can be nicely illustrated 
by the observation that new insights into the nature of 
the weak interaction usually were closely followed by 
the recognition of their importance in some astrophys- 
ical context. Shortly after Pauli postulated the exis- 
tence of the neutrino and Fe r mi de veloped the first theory 
of weak interaction (Fermi, 1934), Gamow and Schoen- 
berg speculated about the possible role of neutrinos in 
stellar evolution and proposed their production in the 
star as an important source for stellar energy losses 
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the bremsstrahlung radiation of neutrino pairs by elec- 
trons would be a very ef fective stellar energy loss mecha- 
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nism ( Pontecorvq 195S). Just afte r the disco v ery o f neu- 
tral weak c urrent |Frecdmari ( 1974 ), Mazurek ( 1975 ), and 
Sate ( 1975] ) recognized that this interaction would result 
in a sizable elastic scattering cross section between neu- 
trinos and nucleons, leading to neutrino trapping during 
the core collapse of a massive star in a type II supernova. 
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The unified model of clectroweak interaction ( Glashow 
1970| ; palarrj , |196cft |Weinberg| , |1967| ) allows deriva 
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tion of accurate cross sections for weak processes among 
elementary particles (i.e. electrons, neutrinos, quarks), 
but also for neutron and protons if proper formfactors are 
taken into account which describe the composite nature 
of the nucleons. However, the situation is different for 
weak interaction processes involving nuclei. Clearly, the 
smallness of the weak interaction coupling parameter al- 
lows treatment of these processes in perturbation theory, 
reducing the calculation basically to a nuclear structure 
problem. However, it has been the inability to adequately 
treat the nuclear many-body problem, which has - and 
in many cases still does - introduced a substantial uncer- 
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tainty into some of the key weak interaction rates used in 
astrophysical simulations. However, the recent few years 
have witnessed a tremendous progress in nuclear many- 
body theory, made possible by new approaches and novel 
computer realizations of established models, but also by 
the availability of large computational capabilities. This 
progress allowed calculation of the rates for many of the 
stellar weak-interaction processes involving nuclei with 
significantly improved accuracy or for the first time. To 
actually know that the calculations are more reliable, im- 
plies the availability of experimental data which test, con- 
strain and guide the theoretical models. Thus, the ad- 
vances in modelling nuclear weak-interaction processes 
in stars also reflects the progress made by experimen- 
talists in recent years which have succeeded to measure 
data which are relevant for the astrophysical applications 
discussed in this review either directly, e.g. half-l ives for 
some short- lived nuclei on the r-process path (Pfciffcr 
et a0, |2001 ) , or indirectly like the Gam ow- Teller distri- 
butions for nuclei in the iron mass range (Osterfeld, 1992) 
which decisively constrain the nuclear models. Another 
recent experimental first has been the measurement of 
charged- and neutral-current neutrino-nucleus cross sec- 
tions. 



This review will report about progress in modelling 
nuclear weak-interaction processes and their possible im- 
plications for stellar evolution and nucleosynthesis. We 
will restrict ourselves to advances achieved by improved 
nuclear models, treating the weak interaction within the 
standa rd model. O f cour se, it has long been recognized 



(e.g. pato and Sato| , |1975| ) that sta rs can be use d as lab- 
oratories for fundamental physics ( jRaffclt , 1996) search- 
ing for new weakly interacting particles or constraining 
exotic components of the weak interaction outside the 
standard model. This field is rapidly growing (see f or ex 



amp le porsico et al\, 2001; Dommguez et 
fclttl |1996| , |1999j . |200q ) . 



199S ; Raf- 



Our review is structured as follows. Following a very 
brief discussion of the required ingredients of the weak 
interaction we introduce the nuclear many-body mod- 
els which have been used in the studies of the weak- 
interaction processes (section |0|). The remaining sections 
are devoted to the results of these calculations and their 
applications to astrophysics which include the solar nu- 
clear reaction network and neutrino problem, the core 
collapse of massive stars, s- and r-process nucleosynthe- 
sis, neutrino nucleosynthesis, explosive hydrogen burning 
and type la supernovae. 

Although generally quite important, weak-interaction 
processes constitute only a part of the many nuclear re- 
actions occuring in stars. For recent reviews about other 
stellar nuclear reactions networks and nucleosynthesis 
the reader is referred to the comprehensive and compe- 



II. THEORETICAL DESCRIPTION 
A. Weak interactions in nuclei 
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FIG. 1 Semileptonic weak processes that occur during the 
evolution of stars. For each process the hadronic current is 
on the left and the leptonic current to the right. The dashed 
circle indicates a bound electron in the initial or final state. 
The four-momentum transfer q\ — (—u>,q) for each pro- 
cess is given in terms of the charged lepton four-momentum 
k\ = (e, k) and the neutrino four-momentum v\ — (y,v)- 
uj, e, and v represent the energy transfer, lepton energy and 
neutrino energy, respectively. In the case of antiparticles the 
directions of the momenta are show as an arrow close to the 
four-momentum label. The first row shows the usual decay 
modes in the laboratory. The second and third rows show 
processes that can occur under stellar conditions. 

Processes mediated by the weak interaction in stars 
can be classified as leptonic (all interacting particles 
are leptons) and semileptonic (leptons interact with 
hadrons via the weak interaction). Leptonic processes 
can be straightfor wardly computed using the standard 
electroweak model (Grotz and Klapdoi, 1990). The cal- 
culation of semileptonic processes (i.e. neutrino-nucleus 
reactions, charged-lepton capture, and /3-decay) is more 
complicated due to the description of the nuclear states 
involved. Fortunately the momenta of the particles turn 
out to be small compared with the masses of the Z, W 
bosons. Thus it is sufficient to consider the semileptonic 
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processes of interest in the lowest-order approximation 
in the weak interaction. Then the interaction can be de- 
scribed by a current-current Hamiltonian density: 



H(x) 



G 



(1) 



where G — GfV uc i for charge-current processes and 
G = Gf for neutral current processes, with Gf the 
Fermi coupling constant and V u d the up- down entry o f 



the C abibbo-Kobayashi-Maskawa matrix ( Groom et 
2000| ). and J^(x) are the weak lcptoni c and 



hadronic density operators (Donnelly and Pcccei, 1979 



Walecka 



1975 



1995). The structure of the lcptonic cur- 



rent j^ix) for a particular process is given by the stan- 



dard electrowea k mod el (Glashow et al, 197C; 3alam 



1968; Weinberg , 1967), and contains both, vector and 



axial-vector components. The standard model describes 
the hadronic current in terms of quark degrees of free- 
dom. Since we are only interested in the matrix el- 
ements of Jfi(x) in nuclei we need only to retain the 
pieces which involve u and d quarks. (The contribution 
from strange quarks is no rmally neglected, but see the 
discussion in section V.C.) As in nuclear physics the 
nucleons are treated as elementary spin-1/2 fermions, 
the Standard Model current is not immediately appli- 
cable. Moreover, nucleons in nuclei interact also via the 
strong interaction. It is then convenient to define an ef- 
fective hadronic current using arguments of Lorentz co- 
variance and isospin invariance of the strong interaction. 
The effective hadronic current can be decomposed into 
strong isoscalar (T — 0) and isovector (T = 1) compo- 
nents and contains both vector (V) and axial-vector (A) 
pieces. The weak charge-changing current is isovector 
with Mt = ±1 and can be written in a general form as: 



J» = V L 



1M T 



A 



lMj 



(2) 



This current governs processes (see figure [T]) such as 
/J^-decay, e-capture, neutrino (z//,Z~) and anti-neutrino 
(j>z, l + ) reactions {I — e, /x or r) . Under the conserved vec- 
tor c urrent (CVC) hypothesis ( Feynman and Gcll-Mann , 
1958| ) the current V^ Mt has a structure identical to the 



isovector part of the electromagnetic current. As a conse- 



quence of this hypothesis the weak charge-changing vec- 
tor current is a conserved quantity. For the weak neutral 
current one has Mt = and, in general, both T = and 
T = 1 pieces can occur. The general form of this current 
is 



J, = ffiv? + Pa'A™ + Pv%» + ft'A™. (3) 
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Assuming that the coupling constants are given by the 
Standard Model we have: 0$ = -2sin 2 9 W , ( T = 0, 

.,2 



1-2 sin" 6 W , { a 



1 ( Donnelly and Peccei , 1979 



p. 25). 0w is the weak mixing angle. The neutral current 
describes weak interactions such as neutrino (y, v') and 
anti-neutrino (D, v') scattering. 

The nuclear transitions that are induced by such weak 
currents (operators) involve initial and final states that 
are usually assumed to be eigenstates of angular momen- 
tum, parity, as well as isospin. It is then convenient to 
do a multipole expansion of the current operators. In 
that way one obtains the Coulomb, longitudinal, trans- 
ver se electric and t ransverse magnetic multipoles defined 
in ( Walecka , 1975 , p. 136). The expressions necessary 
for the calculation of the pro cesses shown in figure l] can 
be obtained from references (Donnelly and Peccei, |1979 



Walecka, 1975| ) in terms of the multipole operators. In 



general the multipole operators are ^4-body nuclear op- 
erators (with A the nucleon number). In practice, at the 
energy scales we are interested in, weak interactions per- 
turb the nucleus only slightly, so that to a good approxi- 
mation one-body components dominate most of the tran- 
sitions. Two-body meson exchange cu rrents and other 
many body effects are negle c ted (see 



Marcucci et al. 



200 1| ; pchiavilla and Wiringaj , [2002| , for a description of 
the nuclear current including two-body operators). It is 
further assumed that a nucleon in a nucleus undergoing a 
weak interaction can be treated as a free nucleon, which 
for the purpose of constructing interaction operators sat- 
isfies the Dirac equation. This latter approximation is 
known as the impulse approximation. For a single free 
nucleon, we have, using Lorentz covariance, conservation 
of parity, time-reversal invariance, and isospin invariance, 
the following general form for the vector and axial vector 
currents 



(k'\';l/2m t ,\Vj' MT \k\-A/2m t ) =iu[k'X') [FP^ + F^'a^ u(k\)(l/2m t ,\I* lT \l/2m t ), 
(fc'Y; l/2m t ,\A™ T \k\; l/2m t ) = iu(k'X') [^757^ - ^FW] u{k\)(l/2m t ,\4 lT \l/2m t ). 

I 
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Here, the plane-wave single nucleon states are labelled transfer, = q 2 — uj 2 , with q — \q\, is defined in fig- 
with the three- momenta fc (A;'), helicities A (A'), isospin urc [l]. Bold letters denote the three- momentum. The 
1/2 and isospin projections m t (m t <). The momentum sing le-nucleon form factors F { P = F ( p(q 2 A, T = 0, 1, 
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X = 



axial, and 



pscudoscalar) are all functions of q 2 u ( 


Beise and McKc- 


own, 


1991; Donnelly and PecccJ, 


1979; 


Kuramoto et al, 


199C; Musolf and Donnelly, 1992) 


. Second class currents 



are not included in equation (j4j) . 
in equations (Q) is contained in 



The isospin dependence 



I*f T = -x 



t±i = T-j^(n ±it 2 ) 



T = 0, M T = 
T = 1,M T = 
T = 1,M T = ±1 



(5) 



To evaluate weak-interaction processes in nuclei, one 
needs matrix elements of the multipole operators be- 
tween nuclear many-body states, labeled | JiMj i ; TjM^) 
which are complicated nuclear configurations of protons 
and neutrons. Using the Wigner-Eckart theorem we can 
write the matrix el ement of an arbi trary multipole oper- 
ator Tjmj-tMt as (Edmonds, 196C) 



(JiMj 1 ;T 1 M Tl \f JMj . TMT (q)\J 2 Mj 2 ;T 2 M T2 ) = (-l)^-^i 

I 



Ji J J 2 

-M, h Mj M, h 

Tt T T 2 

-M Tl M T M T2 



(Jii^lllfjjT^lllJasTa) (6) 



where the symbol 1 1 1 denotes that the matrix element is 
reduced in both angular momentum and isospin. If we 



J 



assume that the multipole operators are one-body oper- 
ators, we can write (Heyde, 1994) 



(Ji;Ti|||f J; r(g)|||J 2 ;T a )=5^ 



a>u]j;T\\\J 2 '-,T't 



V /(2J+1)(2T+1) 
I 



21 , 

a 



\TM*)\\W) 



(7) 



with the sums extending over complete sets of 
single-particle wavefunctions a = n,l,j. The 
tensor product involves the single-particle cre- 



ation operator c 

(_l)Ja-Wa(_l)V2- 



= a a;m ja n Ha and ^ Q = 

a a -- mja -m ta -, witn a the 
destruction operator. The phase factor is introduced 
s o that the oper ator a transforms as a spherical tensor 
( [Edmonds| , [L960D . 

In practice the infinite sums in equation (Q) are approx- 
imated to include a finite number of (hopefully) domi- 
nant terms. The number of terms to include depends 
both of the computed observable and the model used 
(Shell-Model, Random phase approximation, . . . ) . Typ- 
ical nuclear models are non-relativistic, requiring a non- 
relativistic reduction of the single-particle operators; the 
r espec tive e xpressions arc given f o r exam p le by Walecka 
( 1975 ) a nd Donnelly and Pcccci ( 1979 ). Donnelly and 
Haxton (1979) give the expressions for the single-particle 
matrix elements of these operators with harmonic oscilla- 



tor wave functions. Donnelly and Haxtor (198C) provide 



expressions for general wave functions. 

The above discussion presents the general theory of 
semileptonic processes. However, in many applications 
the momentum transfers involved are small compared 



with the typical nuclear momentum Q w -R -1 , with 
R the nuclear radius. In that case, the above formu- 
las can be expanded in powers of (qR) (long-wavelength 
limit) and one obtains the standard approximations to 
allowe d (Gamow- Teller and Fermi) and forbidden transi- 



tions (Bambynek et al, 1977; Behrcns and Biihring, 1971 



1982). In these limits the effect of the electromagnetic 
interaction on the initial or final charged lepton, that has 
been negle c ted in the above expressions, can be included 



( ^chopper , 196(f ) 



B. Nuclear models 

As discussed in the previous section one of the basic 
ingredients for the evaluation of weak-interaction pro- 
cesses involving nuclei is the description of the nuclear 
many-body states. Moreover, the calculation of weak 
processes in stars have to account for the peculiarities of 
the medium (high temperatures and densities) and the 
presence of an electron plasma. When the temperatures 
and densities are small (for example during the r- and 
s-processes) weak transitions could be determined us- 
ing the experimentally measured half-lives (in some cases 
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one has to account for the presence of low lying isomeric 
states). However, as many of the very neutron-rich nu- 
clei that participate in the r-p rocess, are no t curre ntly 



accessible experimentally (see Pfeiffer et ai, 2001, for 



recent experimental advances in the study of r-process 
nuclei), the necessary nuclear properties have to be ex- 
tracted from theoretical models. 

As the degrees of freedom increase drastically with 
the number of nucleons, models of different sophistica- 
tion have to be chosen for the various regions in the nu- 
clear charts. Exact calculations using realistic nucleon- 
nucleon interactions, e.g. by Green's Function Monte 
Carlo tech niques, are restricted to li ght n u clei wit h mass 



A < 10 flCaiison and Schiavilla| , |1998| ; gieperj , |2002 



Wiringa et ai, 200C| ). As an alternative, me t hods based 



effective fi^lH thmcy ( jBgagg pt at 



2om 



van Kolrk 



199ilj ) have re cently been m7\. h. „ ,j for v e rvTiplit rm 



SMMC method available for basically unrestricted model 
spaces. While the strength of the SMMC method is the 
study of nuclear properties at finite temperature, it does 
not allow for detailed nuclear spectroscopy. 

The evaluation of nuclear matrix elements for the 
Fermi operator is straightforward. The Gamow- Teller 
operator connects Slater determinants within a model 
space spanned by a single harmonic oscillator shell (Ofku 
space). The shell model is then the method of choice to 
calculate the nuclear states involved in weak-interaction 
processes dominated by allowed transitions as complete 
or sufficiently converged truncated calculations are nowa- 
days possible for such Ohcu model spaces. The practical 
calculation of the Gamow- Teller distribution is achieved 
by adopting the Lanczo s meth od ( Wilkinsor, [1965 ) as 
proposed by 



Poves, 2000 



Whitehcac| ( 198C); (sec also |Langankc and 



Poves and Nowack, 2001) 



2000, 2001: Park et 



clei (A < 3) ( Marcucci et 

2001a| ,|b|) . For heavier nuclei different approximations are 



required. In particular, restricted model spaces are used 
s o that effective interactions and operators are necessary 
(Hjorth- Jensen et al, 1995|) . For medium-mass nuclei 



(A jt 70) the shell model is the method of choice ( Talmi, et al, 2001; Navratil et 



1993] JI This model explicitly treats all two-body corre 
lations among a set of valance particles by a residual 
interaction. By diagonalizing the respective Hamilto- 
nian matrix in the model space spanned by the inde- 
pendent particle states of the valence particles a quite 
satisfactory description of the ground state, the spec- 
trum at moderate excitation energies and the electro- 
magnetic and weak transitions among these states are 



The calculation of forbidden transitions, however, in- 
volves nuclear transitions between different harmonic os- 
cillator shells and thus requires multi-huj model spaces. 
These are currently only feasible for light nucl ei where 
ab initio shell model calculations are possible (Caurier 

, 2000). Such multi-fiw calcula- 



tions have been used for the calculation of neutrino scat- 



tering from 12 C (Hayes and Townei, 2000; Volpe et al 



obtai ned (Caurier et al, 1994; Marti'ncz-Pinedo et 



1997). In recent times due to progress both in computer 



technology and programming techniques shell-model cal- 
culations are now possible in model spaces which seemed 
impossible only a few years ago; i.e. the diagonalization 
codes ANTOINE or NATHAN developed by Etienne Cau- 
rier allow for complete calculations in the p/-shell where 



the maximum dimension currently attained is 2.3 x 
M = slater determinants for a complete diagonaliza- 



2000). However, for heavier nuclei one has to rely on 
more strongly truncated nuclear models. As the kine- 
matics of stellar weak-interaction processes are often such 
that forbidden transitions are dominated by the collec- 
tive resp onse of the n ucleus the Random Phase Appox- 
imation ( Rowe , 196S ) is usually the method of choice 
(figure ^). Another advantage of this method is that, 
in contrast to the shell model, it allows for global cal- 
culations of these processes for the many nuclei often 
involved in nuclear networks. An illustrative example is 
the evaluation of nuclear half-lives based on the calcu- 
lation of the GT st rength function within t he Q u asipar- 
ticle RPA model (|Krumlinde and Moller| , |l984| ; |Mollci 



and Randrup, 1990). The RPA method considers the 



tion of 60 Zn ( |Caurier| , [2002| ; jMazzocchi et al\ , $00l\ ) 



To treat even larger model spaces in the diagonaliza- 
tion shell model, different schemes are required. One 
method is to expand the nuclear many-body wave func- 
tions in terms of a few symmetry-projected Hartree- 
Fock-Bogoliubov (HFB) type quasiparticle determinants 
( |5chmid| , |2001[ [Bchmid et al\, |1987| |1989|) . A novel ap- 
pro ach, int roduced by [Honma et~al\ ( |T995| ) (see also |5t 



residual correlations among nucleons via one particle one 
hole (lp-lh) excitations in large mnlti-huj model spaces. 
Compared to the shell model, the neglect of higher-order 
correlations renders the RPA method inferior for ma- 
trix elements between individual, non-collective states. 
A prominent example is the GT transition from the 12 C 
ground state to the T — 1 triad in the A = 12 nuclei 
(e.g. Engcl et al, 1996). While the shell model is able to 



suk aj et cdi 2001 ), employs stochastical methods to de- 
termine the most important Slater determinants in the 
chosen model space. As an alternative to the diagonal- 



i zation method the Shell Model Monte C arlo (SMMC) about a factor of 2 (Engcl et al, 1996; Kolbe et 



( Johnson et al. , 1992b ; Koonin et al. , 1997 ) allows calcu- 



reproduce the GT matrix element between these state s 
( |Cohcn and Kurath . 1965 ; Warburton and Brown , 1992 ), 
RPA calculations miss an important part of the nucleon 
correlations and o verestimate t hese matrix elements b y 

, [19941), 



lation of nuclear properties as thermal averages, employ- 
ing the Hubbard-Stratonovich transformation to rewrite 
the two-body parts of the residual interaction by inte- 
grals over fluctuating auxiliary fields. The integrations 
are performed by Monte Carlo techniques, making the 



Recent developments have extended the RPA method to 
include the c omplete set o f 2p-2 h excitations in a given 
model space ( Drozdz et ai , 1990 ). Such 2p-2h RPA mod- 
els have, however, not yet been applied to semileptonic 
weak processes in stars. Moreover, the RPA allows for 
the proper treatment of the momentum-dependence in 



RPA 



SM 



FIG. 2 The most commonly used nuclear models for the calculation of weak processes in stars are the Random Phase Aproxi- 
mation (RPA) and the Shell- Model (SM). In the RPA the basis states are characterized by particle-hole excitations around a 
given configuration (typically a closed-shell nucleus) . In the SM, all the possible two-body correlations in a given valence space 
are considered. Excitations from the core or outside the model space are neglected, but his effect can be included perturbatively 
using effective interactions and operators. 



the different multipole operators, as it can be important 
in certain stellar neutrino-nucleus proces ses (see below) , 



and for the inclusion of the continuum (Buballa et 



1991). Detailed studies indicate that standard and con- 



tinuum RPA calculations yield near ly the same result s 



for total semileptonic cross sections ( Kolbe et al , [2000 ) 



This is related to the fact that both RPA versions obey 
the same sumrules. The RPA has also been extended to 
deal with partial occupation of the orbits so that configu- 
ration mixing in the same shell is included schematically 



( |Ko|bc et al\, |1999b|; |Rowe|, |196S|) 



III. HYDROGEN BURNING AND SOLAR NEUTRINOS 

The tale of the solar neutrinos and their 'famous' prob- 
lem took an exciting twist from its original goal of mea- 
suring the central temperature of the Sun to providing 
convincing evidence for neutrino oscillations, thus open- 
ing the door to physics beyond the standard model of the 
weak interaction. In 1946, Pontecorvo suggested ( Pon- 



teco rvo 



Alvarez 



194£, 1991| ) (later independently proposed by 



1949| ) that chlorine would be a good detector 



material for neutrinos and subsequently in the 1950's 
Davis built a radiochemical neutrino detector which ob- 
served reactor neutrinos via the 37 Cl(i' e , e~) 37 Ar reaction 
(Davi£, p955| ). After the 3 He(a, 7) 7 Be cross section at low 
energies had been found to be significantly larger than 
expected (Holmgren and Johnston, 195S ) and, slightly 



later, the 7 Be(jp, 7) 8 B cross sectio n at low energies had 
been measured ( Kavanagh , I960 ), it became clear that 
the Sun should also operate by what are now known as 
the ppll and ppIII chains and in that way generate neu- 
trinos with energi es high enough t o be de t ectab le by a 
chlorine detector (Cameron, 1958; Fowler, 1958). This 
idea was then seriously pursued by Davis, in close col- 
laboration with Bahcall. The observed solar neutrino 
flux turned out to be lower than predicted by the so- 
lar models (the original solar neutrino problem) (Bahcall 



et al. 



1968 



Davis et 



1968), triggering the develop- 



ment of further solar neutrino detectors, initiating the 
field of neutrino oscillation experiments and, after preci- 
sion helioseismology data (Christensen-Dalsgaard. |200S ) 
boosted the confidence in the solar models, finally cul- 
minating in the conclusive evidence for neutrino oscilla- 
tions in the solar flux. A detailed recent review of the 
solar hydrogen burning and neutrino problem is given in 
( |Kirstcn| , |l999j ). 



The Sun generates its energy from nuclear fusion re- 
actions in the pp chain (see table ||), with a small con- 
tribution by the CNO cycle. Several of these reactions 
are mediated by the weak interaction, and hence create 
(electron) neutrinos which in the standard solar model 
can leave the Sun unhindered. The predicted flux of solar 
neutrinos on the surface of the earth is shown in figure |^. 
These predictions depend on the knowledge of the rele- 
vant nuclear cross sections at solar energies (a few keV) 
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FIG. 3 The energy sp ectrum of neut r inos p redicted by the 
standard solar model (Bahcall et al., 2001). The neutrino 
fluxes from continuum sources (like pp and 8 B) are given m 
the units of number per cm 2 per second per MeV, while the 
line fluxes (e.g. 7 Be) are given in number per cm 2 per sec- 
ond. The ranges of neutrino energies observable in the various 
detectors are indicated by arrows. The uncertainties in the 
various fluxes are given in percent (courtesy of J. N. Bahcall). 



which, wi th the notable excep tion of the 3 He( 3 He,2p) 4 He 
reaction (Bonctti et al, 1999) which has been measured 
directly in the underground laboratory in the Gran Sasso, 
relies on the extrapolation of data taken at higher ener- 
gies ; ^Aj^diesej^^act^ 



olat ons arc quite mild and appear to be under control 



(Adelberger et ai, 1998). The cross section for the ini- 
tial p + p fusion reaction is so low that no data exist 
and the respective solar reaction rate relies completely 
on theoretical modelling. Nevertheless the underlying 
theory is thought to be under control and the uncer- 
taint y in this important rate is estim ated to be about 
1% ( |Kamionkowski and Bahcali . 1994), based on poten- 
tial model calculations, and it is probably even smaller if 



effective fie l d theory is ap plied (Kong and Ravndal, 2001 
Park et ai . 1998 , 2001a ). In the solar plasma the reac- 
t ion rates are slight l y enhanced due to screening effects 
flDzitko et al[ jggj |Gruzinov and BahcaUj |199S| ). The 
most significant plasma modification is found for the lifc- 



time of 'Be with respect to electron capture where cap- 
ture of continuum and bound electrons with the relevant 



scre ening correctio n s have to be accounted for (Gruzinov 



and Bahcall , 1998 ; Johnson et al. , 1992a ) . It is gener- 



Reaction 


Term 


v Energy 




\ /0 ) 


(MeV) 


p + p^ 2 H + e+ +u e 


99.96 


< 0.423 


or 






p + e~ +p — > 2 H + u e 


0.44 


1.445 


2 H+p^ 3 H + 7 


100 




3 He + 3 He -» a + 2p 


85 




or 






3 He + 4 He -> 7 Be + 7 


15 






1 5 < 


f 0.863 90% 
[ 0.385 10% 


7 Li + p -» 2a 






or 






7 Be + p -> 8 B + 7 


0.02 




8 B 8 Be* +e+ + v e 




< 15 


8 Be* -> 2a 






or 






3 He + p -> 4 He + e+ + v e 


0.00003 


< 18.8 



TABLE I The solar pp chains. The neutrino termin ations 
are from the BP2000 solar model ( [Bahcall et aZj |200l|). Th e 
neutrino energies include the solar corrections ( Bahcall 1l997 ) . 



recent direct measurement with special emphasis on the 
control and determina tion of the potential e r rors yielded 



a slightly larger value (Junghans and Snovei, 2002; Jung 



hans et ai, 2001). The neutrino energy distribution aris- 
ing from the subsequent 8 B decay has been measured 
precisely by Ortiz et al. (200C). In principle, high-energy 
neutrinos are also produced in the 3 He+p fusion reaction 
which, however, occurs only in a weak branch in the solar 
pp cycles. Although the calculation of this cross section 
represents a severe theoretical challenge, it appears to be 



of-the-art few-body methods ( 


Kong and Ravndal, 


2001 


Marcucci et ai, 


. 


2001; 


Park et al, 


2001a 


w 



The solar nuclear cro ss sections have been reviewed by 
Adelberger et al. ( 1998| ), including also the reactions oc- 
curing in the CNO cycle. Except for some discrep ancies 
in the 14 N(p,7) 15 Q cross section at low energies (Adel- 



berger et all 1998; Angulo and Descouvemont , 2001), all 



ally believed that the 7 Bc(p, 7) 8 B reaction is the least tor (Bahcall, 1989 



knolvn nuclear input m nuclear models. Altnougn this 
reaction occurs in the weak ppIII chain, the decay of S B 
is the source of the high-energy neutrinos observed by 
the solar neutrino detectors. Recent direct and indirect 
experimental methods have improved the knowledge of 
the 7 Be(p, 7) 8 B rate considerably (Davids et al. (2001) 
and references therein). While these data point to an 
astrophysical S-factor in the range of 18-20 eV b, a very 



relevant solar rates are sufficiently well known 

There are currently 5 solar neutrino detectors oper- 
atin g. Thre e of th em, the homestake chlorine detec- 
' ' P 



et al, 1992), and SAGE ( Abdurashitov et al. 



4 87), GALLEX 1 (|Anselmann 



can only observe charge-current (electron) neutrino re 
actions, while the two water Cerenkov detectors (Super- 



The Gallex detecto r has been recently up graded and has changed 
its name to GNO ( Altmann et al\ , 200C ) 
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Kamiokande flFukuda et~al\ , [l998b| ), SNO ( [Boger et al. 
2000| ) ) also observe neutral-current events, which can be 
triggered by all neutrino flavors. All neutrino detectors 
have characteristic energy thresholds for neutrino detec- 
tion, dictated by the various observation schemes; i.e. 
the detectors are blind for neutrinos with energies less 
than the threshold energy E t h- The pioneering chlorine 
experiment of Davis uses the 37 Cl(z/ e , e~) 37 Ar reaction 
as detector, with E t h = 814 keV. Gallex and Sage de- 
tect neutrinos via 71 Ga(^ e , e~) 71 Ge with the threshold 
energy £^=233.2 keV. In Super-Kamiokande (SK) solar 
neutrinos are identified by the observation of relativistic 
electrons produced from inelastic v + e~ scattering. Due 
to high background at low energies, the observational 
threshold is set to ~ 7 MeV. SNO has an inner vessel 
of heavy water, surrounded by normal water. Like SK, 
this detector can also observe neutrinos via inelastic scat- 
tering off electrons. Additionally, and more importantly, 
SNO can also detect neutrinos by the dissociation of the 
deuteron in heavy water, with the threshold energy of 
order 6 MeV. 

The threshold energies and the predicted solar neutrino 
fluxes are shown in figure One notes that SK and 
SNO are only sensitive to ^3 neutrinos (neglecting the 
weak hep flux), the chlorine experiment detects mainly 



B (70% of Lhe predicted flux by Bahr all el nt -jZdffif) 



and 1 7 Be 



-ob= 



neuLiinus, wlnle Uallex and L?age can alsu 

serve neutrinos generated in the main solar energy source, 
the p+p fusion reaction (54% of the predicted flux). It is 
important to note that the solar neutrino detectors have 



beeii calibialed, 



knuwii iieuliiiiu sum ees. 



using 

The original solar neutrino problem constitutes the 
fact that the earthbound detectors observe less neutri- 
nos than predicted by the solar model. The current com- 
parison is depicted in figure ^. Importantly, Sage and 
Gallex, in close agreement to each other, observe at least 



a ntkutrmo tlux which is consistent with the tact that 
the current solar luminosity is powered by the p + p tu- 



sion reaction ( Abdurashitov et al, 1999; Hampel et al 
1999] ). With improved input (nuclear reaction rates 



opacities, etc.) the solar models evolved and, as a mile- 
stone, passed the stringent test of detailed comparison 
to th e soundspced distribution derived f rom helioseismol- 
ogy ( phristensen-Dalsgaard et al , 1996). It became clear 
that the solution to the solar neutrino problem pointed 
to weak-interaction physics beyond the standard model. 
This line of reasoning was sup ported by the observation 
(Hccger and Robertson, 1996| ) that any solar model as- 
suming standard weak-interaction physics leads to con- 
tradictions between the observed fluxes in the various 
detectors. 



It ha s been speculated already for a long time (Pan 
1968] ) that the solution to the deficient observed 



tecoivo 



neutrino flux lies in the possibility that neutrinos change 
their flavor on their way from the center of the Sun to the 
earthbound detectors. Neutrino oscillations can occur if 
the flavor eigenstates (the physical z/ e ,f M ,j/ r neutrinos) 
are not identical with the mass eigenstates (yi, v%, v-^) of 



the weak Hamiltonian, but rather are given by a unitary 
transformation of these states defined by a set of mix- 
ing angles. Importantly, oscillations between two flavor 
states can only occur if at least one of these states does 
not propagate with the speed of light implying that this 
neutrino has a mass different from zero; more precisely 
Am 2 = rrvy — m 2 ^ 0, where 772,1,2 are the masses of the 
oscillating neutrinos. As all neutrino masses are assumed 
to be zero in the Weinberg-Salam model, the observation 
of neutrino oscillations opens the door to new physics be- 
yond the standard model of weak interaction. Neutrino 
oscillations can occur for free-propagating neutrinos (vac- 
uum oscillations). However, their occurence can also be 
influenced by the environment. In particular, it has been 
pointed out that the high-energy (v e ) solar neutrinos can, 
for a certain range of mixing angles and mass differences, 
transform resonantly into other flavors, mediated by the 
interaction of the v e neutrinos with the electrons in the 
solar plasma, resulting in matter-enhanced oscillations 



(the so-calle d MSW effect |Mikheyev and Smirnovj , |1986 
Wolfenstein|, |1978|). 



First clear evidence for neutrino oscillations was re- 
ported by the SK collaboration which observed a deficit 
of ^-induced events from atmospheric neutrinos and 
could link this deficit to v u <^ v T oscillations (Fukuda 



et al, 1998a; Fukuda et al, 1999| ). [Further evidence for 



neutrin o oscillations has been g i ven b y the L SND collabo- 
ration ( Athanassopoulos et al. , 1998 , |1996| ). This result, 
however, was, for most of the allowed para meter space 
not confirmed by the Karmen experiment (Armbruster 
1998a| ,|b|). The complete LSND result will be tested 



et al. 



by the MiniBoone 2 experiment which is currently under 
construction.] A clear link between neutrino oscillations 
and the solar neutrino problem has been presented re 



cently by a combined analysis (Ahmad et al, 2001) of 



the fi r st SN O data with the precise SK data (Fukuda 
et al. , 2001 ). SNO measured the integrated event rate 
above the kinetic energy threshold T G ff = 6.75 MeV (the 
electron energy threshold then is E t h = T c ft + 0.511 
MeV) for charged-current (CC) reactions on the deuteron 
and inelastic electron scattering (ES). As no evidence 
for a deviation of the spectral shape from the predicted 
shape under the no-oscillation hypothesis has been ob- 
served, the integrated rate could be converted into the 
measured 8 B neutrino flux, resulting in $3^0 i Ve ) = 
1.75±0.07(stat)±0.12(syst)±0.05(theor)xl0 6 cm" 2 s~\ 

^SNO 1 



v) = 2.39±0.34(stat)±0.15(syst) x 10 6 cm" 2 s 



The SNO electron scattering flux result agrees with 
the more precise measurement from SK which yields 
$fg(i/) = 2.32 ± 0.03(stat) ± 0.08(syst) x 10 6 cm" 2 s" 1 . 
We note that the charged-current reaction can only be 
triggered by v e neutrinos at the energies of the solar 
neutrinos. Thus, from the measurement of the v e + D 
event rate, SNO has determined the solar i/ e -flux arriv- 
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FIG. 4 Comparison of the predicted solar neutrino fluxes (Bahcall et al, 2001) for the various neutrino detectors with the 
observed fluxe s. The unit for the chlorine detector (CI) and the Gallium detectors (Ga) are Solar Neutrino Units (SNU) (see 
Bahcall, 1989), while for SNO and SK the comparison is in percentage of the predicted flux (courtesy of John N. Bahcall). 



ing on earth stemming from the decay of 8 B. On the other the same conclusions as the earlier SNO results (Ahmad 



hand, neutrino-electron scattering can occur for all neu- et al, 2001) showing a clear excess of neutral-current 



trino types, whereby the v e + e~ cross section is about 
seven times larger than the z^ jT + e~ cross section. If no 
oscillations involving solar v e neutrinos occur, the SNO 
charged-current flux ^g^o an< ^ the SK inelastic electron 
scattering flux should be the same; that is excluded 
by 3.3 a. [The exclusion is even slightly more severe if 
the recent revision of the v e + D cross section includ- 



ing radiative corrections is considered (Kurylov et al 
2002j ) .] If v e v^ T oscillations occur, <I>g£ should be 
larger than ^g^o as ^ then contains additional neutral- 
current contributions from v^ r neutrinos. From a best 
fit to the SNO and SK data (see figure ||), this con- 
tribution has been determined as ( with la uncertainty ) 
$„ T = 3.69 ± 1.13 x 10 6 cm" 2 s" 1 ( [Ahmad et al\ |2001[ ), 



implying that the total solar flux is (i/ e ) + 3>(v lliT ) 
5.44±0.99 x 10 6 cm -2 s _1 . This result agrees very nicely 
with the 8 B neutrino flux predicted by the solar model 



(Bahcall et al, 2001) (5.05 x 10 b cm~ 2 s" 1 ). 



For years the measurement of the neutral-current v+D 
reaction at SNO has been anticipated as the 'smoking 
gun' for solar neutrino oscillations. After finishing this 
review, the first results of this milestone experiment have 
been published (Ahmad et al, 2002a). They lead to 



over charged-current events, as expected if neutrino os- 
cillations are the origin of the solar neutrino problem. 
Furthermore, the observed neutral-current event rate is 
again consistent wit h the prediction of the solar model 



(Bahcall et al, 2001) 



A global analysis of the latest solar neutrino data 
including the SNO charged-current rate favors matter- 
enhanced neutrino osci ll ation s with large mixing angles 
( Krastev and Smirnov| , [2002 ). Considering the recent 
constraints on the 7 Be(p, 7) 8 B cross section and the re- 
spectively predicted 8 B solar neutrino flux, vacuum os- 
cillations are essentially ex c luded . A similar result is 
obtained by Bahcall et al. (2002) including the recent 



day-ni ght asymmetry measured at SNO (Ahmad et 
2002b[ ) 



For more than 30 years the solar neutrino problem has 
been a demanding challenge for experimentalists and the- 
orists, for nuclear, particle and astrophysicists alike. The 
challenge appears to be mastered, leading to new physics 
and without the need of the many desperate solution at- 
tempts put forward over the years. 
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FIG. 5 Flux of v^, T neutrinos vs. electron neutrinos as de- 
duced from the SNO and SK 8 B neutrino data. The diagonal 
bands show the total 8 B flux as predicted by the sta ndard so- 
lar model (SSM, dashed lines) (Bahcall et al., 2001) and that 
derived from the SNO and SK measurements (solid li nes) 



The interc e pts w ith the axis represent la errors (from Ah- 

" pool. 



mad et al. 



IV. LATE-STAGE STELLAR EVOLUTION 

A. General remarks 

Weak interactions play an essential role already during 
hydrostatic burning. Its importance lies in the fact that 
the neutrinos generated by these processes can leave the 
star unhindered, thus carrying away energy and hence 
cooling the star. While the consideration of energy losses 
by neutrinos is already required during hydrogen burn- 
ing (see above) , the heat flux in the early stages of stel- 
lar burning is predominantly by radiation. This changes, 
following helium burning, when the stellar temperatures 
reach ~ 5 x 10 s K and neutrino-antineutrino pair produc- 
tion and emission becomes the leading energy loss mech- 
anism. The respective cooling rate is a local property of 
the star depending on density p and, very sensitively, on 
temperature T; i.e. the energy loss rate for vv emission 
scales approximately like T , implying that the hot in- 
ner regions of the star cool most effectively. However, 
the dominant nuclear reactions, occuring after helium 
burning, have even stronger temperature dependences. 
For example, the heat production e in the 12 C+ 12 C or 
16 0+ 16 fusion reactions, which dominate hydrostatic 
carbon and oxygen burning, scales like e w T 22 and w T 35 
around T = 10 9 K. As a consequence of the tempera- 
ture gradient in the stellar interior and the vast difference 
in the temperature sensitivity, nuclear reaction heating 
overcomes the neutrino energy loss in the center. How- 
ever, in the cooler mantle region surrounding the core 
neutrino cooling dominates. The resulting entropy diffcr - 
ence leads to convective instabilities (see Arnett, 1996 ). 
First attempts of modelling late-stage stellar burning and 
nucleosynthesis including a two-dimensional treatment of 



The importance of convection has, of course, al- 
ready been noticed before and is accounted for in one- 
dimensi onal models w ithin the so-called mixing-length 
theory ( Clayton , 1968 ). It has been found that this con- 
vective transport is far more efficient at carrying energy 
and mixing the matter composition than radiation trans- 
port. For example, convection dominates the envelope 
region in massive stars during helium shell burning, as 
can be seen in fi gure ^| which shows t he en ergy history 
of a 22 M Q star flHcger and Woosleyj , |200l| ). The figure 
also identifies the various subsequent energy reservoirs of 
the star: hydrogen, helium, carbon, neon, oxygen, and 
silicon core and shell burning. However, the figure also 
demonstrates the importance of neutrino losses which, 
following oxygen core burning, can overcome the nuclear 
energy generation, except at the high temperatures in the 
very inner core. Obviously weak-interaction processes are 
crucial in this late epoch of massive stars. This is not only 
true for the star's energy budget, but these processes can 
also alter the matter composition and entropy which, in 
turn, can affect the location and extension of convective 
shells, e.g. during oxygen and silicon burning, with subse- 
quent changes in the stellar structure. Such effects have 
recently been observed after the improved shell-model 
weak-interaction rates (subsection III. 2) have been in- 
corporated into stellar models, (figure || is already based 
on these rates.) 



B. Shell-model electron capture and 8 decay rates 

The late evolution stages of massive stars are strongly 
influenced by weak interactions which act to determine 
the core entropy and electron to baryon ratio, Y e , of the 
presupernova star, hence its Chandrasekhar mass which 
is proportional to Y 2 . Electron capture reduces the num- 
ber of electrons available for pressure support, while beta- 
decay acts in the opposite direction. Both processes gen- 
erate neutrinos which, for densities p < 10 11 g cm~ 3 , 
escape the star carrying away energy and entropy from 
the core. 

Electron capture and beta decay during the final evo- 
lution of a massive star are dominated by Fermi and 
Gamow- Teller (GT) transitions. While the treatment 
of Fermi transitions (important only in beta decays) is 
straightforward, a correct description of the GT tran- 
sitions is a difficult problem in nuclear structure. In 
the astrophysical environment nuclei are fully ionized so 
one has continuum electron capture from the degener- 
ate electron plasma. The energies of the electrons are 
high enough to induce transitions to the Gamow- Teller 
resonance. Shortly aft e r the discovery of this collective 
excitation 



Bcthc et al. (1979) recognized its importance 



convection is reported in (Baleisis and Arnett, 2001) 



for stellar electron capture. The presence of the degener- 
ate electron gas blocks the phase space for the produced 
electron in beta decay. Then the decay rate of a given nu- 
clear state is greatly reduced or even completely blocked 
at high densities. However, due to the finite temperature 
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FIG. 6 Energy history of a 22 Mq star as a function of time until core collapse. The y-axis defines the included mass from the 
center. Hydrogen and helium core and shell burning are major energy sources. In the later burning stages, following oxygen 
core burning, neutrino losses related to weak processes in the stellar interior become increasingly important and can dominate 
over the nuclear energy production. Convection plays an important role in the e nvelope outside the helium burning shell, but 
also in shells during oxygen and silicon burning (from Heger and Woosley, 200i). 



excited states in the decaying nucleus can be thermally 
populated. Some of these states are connected by strong 
GT transitions to low-lying states in the daughter nu- 
cleus that with increased phase space can significantly 
contribute to the stellar beta decay rates. The impor- 
tance of these states in the parent nucleus for beta-decay 
was first recognized by Fuller, Fowler and Newman (com- 
monly abbreviated as FFN) who coined the term "back- 
resonances" (see figure 0) . 

Over the years, many calculations of weak interaction 
r ates for astrophysical app l ications have become available 
flAufdcrhcide et al[ |l994cj; |Hanscn|, |l9G6|, |l9Ggj: [Vlazurek 



mation in the mid 1980s for individual transitions be- 
tween ground states and low-lying excited states in the 
nuclei of interest. Recognizing that this only saturated a 
small part of the Gamow- Teller distribution, they added 
the collective strength via a single-state representation. 
Both, energy position and strength collected in this sin- 
gle state were determined using an independent particle 
model (IPM). 

Recent experimental data on GT distributions in iron 
" H99l ~ 



197$ pVlazurek et al.[ |1974| [Takahashi et al\, |l973|). For et al 



group nucl ei ( Alford et 
198.5|, |199C|; |E1-Kateb et aUTl994|; 



1990; Anderson et 



Rapaport et al. , 1983| 



Ronnqvist et al, 1993; Vetterli et al, 1990; Williams 



approximately 15 years tho ugh, the standard in the field 
has been the ta bulati ons of Fuller, Fowler, and Newman 
( |1980| , |l982aj b] |l985| ). These authors calculated rates 



for electron capture, positron capture, beta-decay, and 
positron emission plus the associated neutrino losses for 
all the astrophysically relevant nuclei ranging in mass 
number from 21 to 60. Their calculations were based 
upon an examination of all available experimental infor- 



_ 1995 ) measured in charge exc hange reactions 
( poodman et al. , 1980| ; Ostcrfeld , 1992 ), show that the 
GT strength is strongly quenched, compared with the 
independent particle model value, and fragmented over 
many states in the daughter nucleus. Both effects are 
caused by the residual interaction among the valence nu- 
cleons and an accurate description of these correlations 
is essential for a reliable evaluation of the stellar weak- 
interaction rates due to the strong phase space energy 



12 



Electron capture 








1.0 






0.8 


on 




0.6 






0.4 


15 




0.2 


, — „, 
> 




0.0 


CD 
















0.6 














Q) 






■*— > 


0.4 




00 




5 


h- 






o 


0.2 






0.0 



(Z+1.A) 

FIG. 7 The figure shows schematically the electron capture 
and beta decay processes in the stellar environment. Elec- 
tron capture proceeds by Gamow- Teller transitions to the 
GT+ resonance. In the case of beta decay both the Fermi 
and Gamow- Teller resonances are typically outside of the Qp 
window, and hence are not populated in laboratory decays. 
Due to the finite temperature in stars excited states in the 
decaying nucleus can be thermally populated. Some of these 
states have strong GT transitions to low-lying states in the 
daughter nucleus. These states in the decaying nucleus are 
called "backresonances" . 



dependence, particularly of the stellar electron-capture 
rates. The shell model is the only k nown tool to reliably 
descr ibe G T distribu tions in nu clei (|Brown and Wildcn ' 
thal[ [1988D . Indeed, paurier et 



(1999) demonstrated 
that the shell model reproduces all measured GT + dis- 
tributions (in this direction a proton is changed into a 
neutron, like in electron capture) for nuclei in the iron 
mass range very well and gives a very reasonable account 
of t ho oxporimontally known GT_ distributions (in this 



1.0 



0.5 - 



0.0 




I ' I ' I 
2 2 4 6 8 10 12 



-51 V 


I _ 


1 V 

r Jf, 


1 ; 


_ 55 Mn l 








iii 
- 59 Co _ 

E 


III 







-2 2 4 6 8 10 12 

E (MeV) 



E (MeV) 

Compariso n of shell model GT+ distribution s with 



FIG. 

experimental data (Afford et al. 



1993 



El-Kateb et al. 



1994 



Ronnqvist et al, 199S) for selected nuclei. The shell model 



results (discrete lines) have been folded with the experimen- 
tal res olution (histograms) . The arrows indicate the positions 
(1982VJ) placed the GT resonance in their 



where Fuller et al 



calcu lations of the stell ar weak-interaction rates (adapted 
from paurier et al. , 199£ ) . 



calculated shell-model rates for all the relevant weak pro- 
cesses for sd-shell nuclei (A — 17-39). This wor k was 
then extended to heavier nuclei (A = 45-65) by Lan 



dire ction, a Beaton is changed into a protoy like in ,6 ganke and Martmcz-Pinedcj ( gOOjJ ) based on shell-model 



decay). Further, the lifetimes of the p/-shell nuclei and 
their spectroscopy at low energies are simultaneously also 
described well. Figure || compares the shell model GT + 
distributions to the pioneering measurement performed 
at TRIUMF. These measurements had a typical energy 
reso lution of ~ 1 MeV. Recently developed techniques. 
invo lving for example ( 3 Hc, t) ( Fujita et al. ,_ 1996 ) and 
(d, 2 He) (Wortche, 2001) charge-exchange reactions at 



calculations in the complete pf shell. Following the spirit 
of FFN, the shell model results have been replaced by ex- 
perimental data (energy positions, transition strengths) 
wherever available. 

Weak interaction rates have also been comp uted using 
the proton-neutron quasipartic lc RPA model ( Nabi and 
Klapdor-Klcingrothaus], |l999a , ^ ) and the spectra l distri - 



bution theory ( |Kar et al , 1994 



Sutaria and Ray , 1995 ) 



intermediate energies, demonstrated in pilot experiments 
an improvement in the energy resolution by an order of 
magnitude or more. Again, the shell model calculations 
agree quite favorably with the improved data. 



Seve ral years ago , |Aufdcrhcidc (1991) and Aufdcrhcidc 
et a| ( |1996| , |l993gj ff 1T994cj ) pointed out that the inter- 



acting shell model is the method of choice for the cal- 
culation of stellar weak-interaction rates. Following the 
work by |Brown and Wildentha] Ql98S| ), pda et al\ ( |1994| ) 



After oxygen burning, the important weak processes 
are electron captures and beta decays on nuclei in the 
iron mass range (A ~ 45-65). Conventional stellar mod- 
els descri bed these weak processes usin g the ra tes esti- 
mated by Fuller, Fowler, and Ncwman| ( 1982b| ). These 
rates are compared to the shell model electron capture 
rates in figure ^ at relevant temperatures and densities. 
Importantly the shell model rates are nearly always lower 
than the FFN rates. Thus this difference represents a sys- 
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FIG. 9 Shell model electron-capture rates as a function of 
temperature (Tg measures the temperature in 10 9 K) and for 
selected densities (p7 defines the density in 10 7 g cm -3 ) and 
nuclei. For comparison, the FFN rates are given by the full 
points. 



tematic trend, which is not expected to be washed out if 
the many nuclei in the stellar composition are considered. 
The difference is caused, for example, by the reduction of 
the Gamow- Teller strength (quenching) compared to the 
IPM value and a systematic misplacement of the Gamow- 
Tcllcr centroid in nuclei with ce r tain pairing structure 
( Langankc and Martmcz-Pinedo , 2000 ). In some cases, 
experimental dat a, which were not available to FFN, bu t 
r.nn Id be used bjqLanganke and Marti"ne7-PineHn| ( [?001 \ 
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FIG. 10 Shell model beta-decay rates as a function of temper- 
ature (Tg measures the temperature in 10 9 K) and for selected 
densities (p-j defines the density in 10 7 g cm -3 ) and nuclei. 
For c omparison, the FFN ra t es ar e given by the full points 
(from 



Martmez-Pinedo et al 



2000) 



and Woosley (WW) used the FFN rates for electron 



capture and an older set of beta decay rates ( Mazurek , 

As a side-remark we note 



1973: Mazurek et 



1974| ) 



that late-stage evolution of massive stars is quite sensi- 
tive to the still not sufficiently well known 12 C(a, 7) 16 
rate . The value adopted in the standard WW and in 
the |Hcger et al.\ models [S(E = 300 keV) = 170 keV b] 
agrees, however, rather nicely with the recent da ta analy- 
sis [5(300) = 165 ±50 keV b flKunz et al\ , |200l| )] and the 
value derived from n ucleosynthesis arguments by Weaver 
led llo significant chafes. The FFN and shell-modelheta and Wooskyj (|993|), [5(300) = 170 ± 20 keV b]. 



Martmez-Pinedc 



decay rates are compared in figure |1^ 
et a|] ( |2000| ) discuss the differences between the two rate 
sets. 



C. Consequences of the shell model rates in stellar models 



illustrates the consequences of the shell 
interaction rates for presupernova models 
three decisive quantities: the central 
the entropy, and the iron core 
at the onset of core col- 



Heger et al. ( 2001a. bj) have investigated the influence 
of the shell model rates on the late-stage evo lution of 
massive stars by r epeating the calculations of Woosley 
and Weaver ( 1995 ) keeping the stellar physics, except for 
the weak rates, as close to the original studies as possible. 
The new calculations have incorporated the shell-model 
weak interaction rates for nuclei wit h mass number s 
A = 45-65, supplemented by rates from Oda et al. (1994) 
for lighter nuclei. The earlier calculations of Weaver 



Figure 
model weak 
in terms of the 
electron-to-baryon ratio Y t 
mass. The central values of Y e 
lapse increased by 0.01-0.015 for the new rates. This is a 
significant effect. For example, a change from Y e = 0.43 
in the WW model for a 20 M star to Y e = 0.445 in the 
new models increases the respective Chandrasekhar mass 
by about 0.075 M & . We note that the new models also 
result in lower core entropies for stars with M < 20 M Q , 
while for M > 20 M Q , the new models actually have a 
slightly larger entropy. The iron core masses are gener- 
ally smaller in the new models where the effect is larger 
for more massive stars (M > 20 M©), while for the most 
common supernovae (M < 20 M©) the reduction is by 
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FIG. 11 Comparison of the center values of Y e (left), the iron core sizes (middle) and the central entropy (right') for 11-40 Mm 
stars between the WW models and the ones using the shell model weak interaction rates (LMP) (from Heger et ai , 2001a ). 
The lower parts define the changes in the 3 quantities between the LMP and WW models. 



about 0.05 Mq. [We define the iron core as the mass inte- 
rior to the point where the composition becomes at least 
50% of iron group elements (A > 48)]. This reduction of 
the iron core mass appears to be counterintuitive at first 
glance with respect to the slower electron capture rates 
in the new models. It is, however, related to changes in 
the entropy profile during silicon shell burning which rc- 
duces the gr o wth of the iron core just prior to collapse 
QHeger et al\ , [2001a| ). 

It is intriguing to speculate what effects these changes 
might have for the subsequent core collapse and super- 
nova explosion. At first we note that in the current su- 



pernova picture (Bcthc, 1990; Burrows, 200C; Langanke 
and Wiescher, 2001; Woosley et ai, 2002) gravitation 



overcomes the resisting electron degeneracy pressure in 
the core, leading to increasing densities and tempera- 
tures. Shortly after neutrino trapping at densities of a 
few 10 11 g cm -3 , an homologous core, which stays in 
sonic communication, forms in the center. Once the core 
reaches densities somewhat in excess of nuclear matter 
density (a few 10 14 g cm -3 ) the nuclear equation of state 
stiffens and a spring-like bounce is created triggering the 
formation of a shock wave at the surface of the homolo- 



gous core ( Bethe , 199C ) . This shock wave tries to traverse 
the rest of the infalling matter in the iron core. However, 
the shock loses its energy by dissociation of the infalling 
matter and by neutrino emission, and it is generally be- 
lieved now that supernovae do not explode promptly due 
to the bounce shoc k. Probably, t his happens in the 'de- 
layed mechanism' (Wilson, 1989) where the shock is re- 
vived by energy deposition from the neutrinos generated 
by the cooling of the proto-neutron star, the remnant in 
the center of the explosion. 

With the larger Y e values, obtained in the calculations 



with the improved weak rates, the core contains more 
electrons whose pressure acts against the collapse. It is 
also expected that the size of the homologous core, which 
scales like ~ Y^ with the Y e value at neutrino trapping, 
should be larger. This, combined with the smaller iron 
cores, yields less material which the shock has to traverse. 
Furthermore, the change in entropy will affect the mass 
fraction of free protons, which in the later stage of the 
collapse contribute significantly to the electron capture. 
For presupernova models with masses M < 20 Mq, how- 
ever, the number fraction of protons is very low (< 10~ 6 , 
2001a| ) so that for these stars electron cap- 



Hcger et 



turc should still be dominated by nuclei, even at densities 
in access of 10 10 g cm -3 . We will return to this problem 
below. 

To understand the origin of these differences it is il- 
lustrative to investigate the role of the weak-interaction 
rates in greater details. The evolution of Y e during the 
collapse phase is plotted in figure |l^. Weak processes 
become particularly important in reducing Y e below 0.5 
after oxygen depletion (~ 10 7 s and 10 6 s before core col- 
lapse for the 15 Mq and 25 Mq stars, respectively) and 
Y e begins a decline which becomes precipitous during sili- 
con burning. Initially electron capture occurs much more 
rapidly than beta decay. As the shell model rates are gen- 
erally smaller than the FFN electron capture rates, the 
initial reduction of Y e is smaller in the new models; the 
temperature in these models is correspondingly larger as 
less energy is radiated away by neutrino emission. 

An important feature of the new models is demon- 
strated in figure 13, Beta decay becomes temporarily 



competitive with electron capture after silicon depletion 
in the core and during silicon shell burning (this had 
been foreseen in Aufderheide et ai, 1994b). The pres- 
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FIG. 13 Comparison of the change of the Y e value with time, 
\dY e /dt\ due to electron capture and beta decay in a 15 Mq 
star. In general, the Y e value decreases with time during the 
collapse, caused by electron captures. The loops indicate that 
during this period beta-decay, which increases Y e , dominates 
over electron capture. 



ence of an important beta decay contribution has two 
effects. Obviously it counteracts the reduction of Y e in 
the core, but equally important, beta decays are an ad- 
ditional neutrino source and thus they add to the cooling 
of the core and a reduction in entropy. This cooling can 
be quite efficient as often the average neutrino energy in 
the involved beta decays is larger than for the competing 
electron captures. As a consequence the new models have 
significantly lower core temperatures than the WW mod- 
els after silicon burning. At later stages of the collapse, 
beta decay becomes unimportant again as an increased 
electron chemical potential drastically reduces the phase 
space. 

We note that the shell model weak interaction rates 
predict the presupernova evolution to proceed along a 
temperature-density- Ye trajectory where the weak pro- 
cesses are dominated by nuclei rather close to stability. 
Thus it will be possible, after next generation radioactive 
ion-beam facilities become operational, to further con- 
strain the shell model calculations by measuring relevant 
GT distributions for unstable nuclei by charge- exchange 
reaction, where we like to point out again that the GT + 
distribution is also crucial for stellar /3-decays. Figure [12| 
identifies those nuclei which dominate (defined by the 
product of abundance times rate) the electron capture 
and beta decay during various stages of the final e volu- 
tion of 15 Mq and 25 Mq stars. Hcgcr et al. ( |2001b ) give 
an exhaustive list of the most important nuclei for both 
electron capture and beta decay during the final stages 
of stellar evolution for stars of different masses. 

In total, the weak interaction processes shift the mat- 
ter composition to smaller Y e values (see Fig. [12]) and 
hence more neutron-rich nuclei, subsequently affecting 
the nucleosynthesis. Its importance for the elemental 
abundance distribution, however, strongly depends on 
the location of the mass cut in the supernova explosion. 
It is currently assumed that the remnant will have a bary- 
onic mass between the iron core and oxygen shell masses 
2002| ). As the reduction of Y e occurs 



(Woosley et al 



mainly during silicon burning, it is essential to determine 
how much of this material will be ejected. Another im- 
portant issue is the possible long-term mixing of material 
during the explosion (e.g. Kifonidis et ai, |2000 ). Changes 



of the elemental abundances due to the improved weak- 
interaction rates are rather small as the differences, com- 
pared to FFN, occur in regions of the star which are 
probably not ejected (however, for type la supernovae, 
see below). The weak interaction also determines the de- 
cay of the newly synthesized nuclei in supernova explo- 
sions. Some of them are proton-rich nuclei that decay by 
orbital electron capture, leaving atomic K-shell electron 
vacancies. The X-rays emitted can escape the supernova 
ejecta for sufficiently long-lived isotopes and can possibly 
be detec t ed by the current generation of X-ray telescopes 
flLeisingt [200 lfr . 

In dense stellar environment the electron capture rates 
have to be corrected for screening effects caused by the 
relativistically degenerate electron liquid. Such studies 
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have be en recently perfo rmed within the linear response 
theory ( Itoh et al. , 2002 ) who find typical screening cor- 
rect ions of order a few percent. 



cient nuclear physics input? Although first steps have 



V. COLLAPSE AND POST-BOUNCE STAGE 

The models, as we have described them above, con- 
stitute the so-called presupernova models. They follow 
the late-stage stellar evolution until core densities just 
below 10 10 g cm -3 a nd temperatures between 5 and 10 



GK. (More precisely, Woosley and Weavci ( 1995 ) define 
the final presupernova models as the time when the col- 
lapse velocity near the edge of the iron core first reached 
1000 km s _1 .) As we have stressed above, stellar evo- 
lution until this time requires the consideration of an 
extensive nuclear network, but is simplified by the fact 
that neutrinos need only be treated as a source for en- 
ergy losses. This is no longer valid at later stages of the 
collapse. As neutrinos will eventually be trapped in the 
collapsing core and their interaction with the surround- 
ing matter is believed to be crucial for the supernova 
explosion, computer simulations of the collapse, bounce 
and explosion necessitate a detailed time- and space- 
dependent bookkeeping of the various neutrino (y e , v^, v T 
neutrinos and their antiparticles) distr ibutions in th e 
core. Built on the pioneering work by Brucnn ( p.985 ), 
this is done by multi-group (neutrinos of di fferent flavor 



and e nergy) Boltzmann neutrino tran sport ( Liebcndorfcr 



et all 2001; Mezzacappa et al., 2001; Rampp and Janka 
200C). Advantageously, the temperature during the col- 
lapse and explosion are high enough that the matter 
composition is given by nuclear statistical equilibrium 
without the need of reaction networks for the strong and 
electromagnetic interaction. The transition from a rather 
complex global nuclear network, involving many neutron, 
proton and a fusion reactions and their inverse, to a 
quasi-statistical equilibrium, in which reactions within 
mini-cycles are fast enough to bring constrained regions 
of the nuclear chart into equilibrium, to finally global 
nuclear statistic al equilibrium is extensively discussed by 
(IWoosleyj p8^) . 

Presupernova models are the input for collapse and 
explosion simulations. Currently, one-dimensional mod- 
els with sophisticated neutrino transport do not ex- 



plode ( 


Liebendorfer et al. 


, 2001 




Mezzacappa et a/., 


2001 


; Rampp and Janka, 


200C), including first attempts 



with the presupernova models derived wit h the improved 



weak - interaction rates discussed above ( Messer et al 
2002). Explosions can, however, be achieved if the 



shock revival in the delayed mechanism is modelled by 



y transfer ( 


Burrows et al. 


Janka and Mullci, 


1996 


)■ 



1995 



Herant et al\ |1994| ; |Janka and Mullerj , ^996|). Thus 
the intriguing question arises: Are supernova explosions 
three-dimensional phenomena requiring convective mo- 
tion and perhaps rotation and magnetic fields? Or do 
one-dimensional models fail due to incorrect or insuffi- 



Janka et al 


3 


2001; 


Woosley 



et al. , |2002j ), these require extremely demanding and 
computationally challenging simulations. In the follow- 
ing we will briefly discuss some nuclear physics ingredi- 
ents in the collapse models and their possible improve- 
ments. 

The c r ucial weak pr o cesses during the collapse ar e 



(Brucnn, 1985; Burrows, 2001; Rampp and Janka, 2002) 
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Here, a nucleus is symbolized by its mass number A 
and charge Z, N denotes either a neutron or a proton 
and v represents any neutrino or antineutrino. In the 
early collapse stage, before trapping, these reactions pro- 
ceed dominantly to the right. We note th at, due to th e 
generally accepted collapse picture (e.g. [Bethe , 1990 ) 



elastic scattering of neutrinos on nuclei (BgJ) is mainly 
responsible for the trapping, as it determines the dif- 
fusion time scale of the outwards streaming neutrinos. 
Shortly after trapping, the neutrinos are thermalized by 
the energy downscattering, experienced mainly in inelas- 
tic scattering off electrons (pR)- The relevant cross sec- 



tions for these processes are readily derived ( Brucnn , 
1985). For elastic neutrino- nucleus scattering one usu- 



ally makes the simplifying assumption that the nucleus 
has a J = + spin/parity assignment, as appropriate for 
the ground state of even-even nuclei. The scattering pro- 
cess is then restricted to the F ermi part of the neutral 
current (pure vector coupling) ( |Freedman , 1974; Tubbs 
and Schramm, [L97|) and gives rise to coherent scattering; 
i.e. the cross section scales with A 2 , except from a correc- 
tion ~ (N — Z)/A arising from the neutron excess. This 
assumption is, in principle, not correct for the ground 
states of odd- A and odd-odd nuclei and for all nuclei at 
finite temperature, as then J > and the cross section 
will also have an axial- vector Gamow- Teller contribution. 
However, the relevant GTo strength is not concentrated 
in one state, but rather fragmented over many nuclear 
levels. Thus, one can expect that the GT contributions 
to the elastic neutrino-nucleus cross sections are in gen- 
eral small enough to be neglected. 

Reactions ( |Sa| ) and d8c| ) are equally important, as they 
control the neutronization of the matter and, in a large 
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portion, also the star's energy losses. Due to their strong 
phase space sensitivity (~ E%), the electron capture cross 
sections increase rapidly during the collapse as the den- 
sity (the electron chemical potential scales like ~ p 1 ^ 3 ) 
and the temperature increase. We already observed 
above that beta-decay is rather unimportant during the 
collapse due to Pauli-blocking of the electron phase space 
in the final state. We also noted how sensitively the elec- 
tron capture rate on nuclei depends on a proper descrip- 
tion of nuclear structure. As we will discuss now, this is 
also expected for this stage of the collapse, although the 
relevant nuclear structure issues are somewhat different. 



A. Electron capture on nuclei 

The new presupernova models indicate that electron 
capture on nuclei will still be important, at least in the 
early stage of the collapse. Although capture on free pro- 
tons, compared to nuclei, is favored by the significantly 
lower Q- value, the number fraction of free protons Y p , i.e. 
the number of free protons divided by the total number 
of nucleons, is qu ite low (Y n ~ 10~ 6 in the 15 M Q presu- 



pernova model of HegeretaZ. , 2001b). This tendency had 



already been observed before, but has been strengthened 
in the new presupernova models, where the Y e values 
are significantly larger and thus the nuclei present in the 
matter composition are less neutron-rich, implying lower 
Q-values for electron capture. Furthermore, the entropy 
is smaller in stars with < 20 Mq, yielding a smaller frac- 
tion of free protons. 

As the entropy is rather low (Bethc et al, 197S| ), most 
of the collapsing matter survives in heavy nuclei. How- 
ever, Y e decreases during the collapse making the matter 
composition more neutron-rich, hence energetically fa- 
voring increasingly heavy nuclei. In computer studies of 
the collapse, the ensemble of heavy nuclei is described 
by one representative which is generally chosen to be the 
most abundant in the nuclear statistical equilibriu m com- 



posi tion. Due to a s imulat ion of the infall phase ( Mezza- 
cap pa and Brucm] , 1993a. b|), such representative nuclei 



are 



YU Zn and 88 Kr at different stages of the collapse ( Mez- 
zaca ppa], |200l| ). 



In current collapse simulations the treatment of elec- 
tron capture on nuclei is schematic and rather simplis- 
tic. The nuclear structure required to derive the capture 
rate is then described solely on the basis of an indepen- 
dent particle model for iron range nuclei, i.e., consider- 



.fk/2 neutrons (Bethc et al. 




1979; 


Bruenn 


. L985: 


Mezza- 


cappa and Bruenn, 1993a, 


cJ). In particular, this model 



predicts that electron capture vanishes for nuclei with 
charge number Z < 40 and neutron number N > 40, ar- 
guing that Gamow- Teller transitions are blocked by the 
Pauli principle, as all possible final neutron orbitals are 
already occu pied in nuclei with N > 40 (closed pf shell) 



and (Z 



36, N — 52), respectiv e ly. It has been pointed 



out ( |Cooperstcin and Wambach , 1984 ) that this picture 
is too simple and that the blocking of the GT transi- 
tions will be overcome by thermal excitation which ei- 
ther moves protons into gg /% orbitals or removes neutrons 
from the pf shell, in both ways reallowing GT transitions. 
In fact, due to this 'thermal unblocking', GT transitions 
again dominate the electron capture on nuclei for tem- 
peratures of order 1.5 MeV (Cooperstein and Wambach, 
1984). An even more important unblocking effect, which 



is already relevant at lower temperatures is, however, ex- 
pected by the residual interaction which will mix the gg « 
(and higher) orbitals with those in the pf shell. 

A consistent calculation of the electron capture rates 
for nuclei with neutron numbers N > 40 and proton 
numbers 20 < Z < 40, including configuration mixing 
and finite temperature, is not yet feasible by direct shell 
model diagonalization due to the large model spaces and 
many states involved. It can, however, be performed in 
a reasonable way adopting a hybrid model: The cap- 
ture rates are calculated within the RPA approach with 
partial occupation formalism, including allowed and for- 
bidden transitions. The partial occupation numbers rep- 
resent an 'average' state of the parent nucleus and de- 
pend on temperature. They are calculated within the 
Shell Model Monte C arlo approach at finite temperature 



(Koonin et al 



interaction. 



1997) and include an appropriate residual 
Exploratory studies, performed for a chain 
of germanium isotopes [Z = 32), confirm that the GT 
transition is not blocked for N > 40 and still dominates 
the elect ron capture proc e ss for s uch nuclei at stellar con- 



ditions ( Langanke et al. , 2001a ). This is demonstrated 
in figure fT4|, which compares electron capture rates for 
78 Ge calculated within the hybrid model with the re- 
sults in the independent particle model (IPM). For this 
nucleus (N — 46) the rate in the IPM is given solely by 
forbidden transitions (mainly induced by 1 _ and 2~ mul- 
tipoles). However, correlations and finite temperature 
unblock the GT transitions in the hybrid model which 
increases the rate significantly. The differences are par- 
ticularly important at lower densities (a few 10 10 g cm~ 3 ) 
where the electron chemical potential does not suffice to 
induce forbidden transitions. 

We note again that many nuclei are present with simi- 
lar mass abundances during the supernova collapse phase 
and that their relative abundances are approximately 
described by nuclear statistical equilibrium. Figure |l4| 
shows the capture rates for several representative nu- 
clei during the collapse phase, identified by the average 
charge and mass number of the matter composition fol- 
lowing the time evolution of a cer t ain ( M — 0.6 Mq) 

200lh . (Forthecon- 



( Fuller, 1982). Such a situation would, for example, occur 
for the two nuclei 70 Zn and 88 Kr with (Z = 30, N = 40) 



mass trajectory ( Liebendorfer et al. 

ditions shown in figure |l4| IJti Kr and v ^Zn are examples of 
representative nuclei at pu = 10.063 and 0.601, respec- 
tively). The general trend of the rates reflects the com- 
petition of the two main energy scales of the capture pro- 
cess: the electron chemical potential /i e , which grows like 
p 1 / 3 during infall, and the reaction Q-value. As the Q- 
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FIG. 14 Electron capture rates on free proton and sev- 
eral representative nuclei as a function of electron chemi 
cal potential fi e (lower panel). The rates hav e been deter 



mine d for conditions during the infall phase ( Liebendorfei , 
200^ ): (T,y e ,pu)=(0.846,0.429,0.102), (1.133,0.410,0.601), 
(1.259,0.400,1.186), (1.617,0.372,4.287), (2.000,0.349,10.063), 
where the temperature T is given in MeV and pn measures 
the density in 10 11 g cm -3 . The electron chemical poten- 
tial increases during the infall, n e ~ p 1 ^ 3 . The middle panel 
compares the capture rates on 78 Ge in the independent par- 
ticle model (IPM), in which Gamow- Teller transitions are 
Pauli-blocked, with the results obtained in the hybrid model 
(SMMC/RPA) which unblocks these transitions due to cor- 
relations and finite temperature effects. The rates have been 
calculated without neutrino final-state Pauli blocking which 
will become important at trapping densities. The upper panel 
shows the time evolution of the number abundances for neu- 
trons Y n , protons Y p , a-particles Y a and heavy nuclei Yh cal- 
culated along the same stellar trajectory for which the rates 
have been calculated at s elected points. These points are 
marked by asterisks (from Sampaio et at, 2002a). 



value is smaller for free protons (Q — 1.29 MeV) than for 
neutronrich nuclei {Q ~ few MeV), the capture rate on 
free protons is larger than for the heavy nuclei. However, 
this difference diminishes with increasing density. This is 
expected because the electron energies involved (for ex- 
ample, the electron chemical potential is fj, e ~ 18 MeV 
at p = 10 11 g cm -3 ) are then significantly higher than 
the Q-values for the capture reactions on the abundant 
nuclei, (i.e. 93 Kr has a Q-value of about 11 MeV). As also 
nuclear structure effects at the relatively high tempera- 
tures involved are rather unimportant, the capture rates 
on the abundant nuclei at the later stage of the collapse 
are rather similar. However, the capture rate is quite 
sensitive to the reaction Q-value for the lower electron 
chemical potentials. To quantify this argument we take 
the point of the stellar trajectory from figure [l4| with the 
lowest electron chemical potential (fi e ~ 8 MeV) as an ex- 
ample. Under these conditions the capture rates on 70 Cu 
and 76 Ga (both nuclei have Q-values around 4 MeV) are 
noticeably larger than for 78 Ge and 72 Zn with Q- values 
around 8 MeV. However, in nuclear statistical equilib- 
rium the relative mass fraction of 72 Zn (about 1.2 x 1CV 2 ) 
is larger than for 70 Cu (4.0 x 1CT 3 ) or 76 Ga (1.8 x 1CT 3 ). 
The most abundant nucleus, 66 Ni, has a mass fraction of 
4.3 x 10~ 2 and a capture rate comparable to Zn. 93 Kr is 
too neutron-rich to have a significant abundance at this 
stage of the collapse. This discussion indicates that the 
most abundant nuclei are not necessarily the nuclei which 
dominate the electron capture in the infall phase. Thus, 
a single-nucleus approximation can be quite inaccurate 
and should be replaced by an ensemble average. 

What matters for the competition of capture on nu- 
clei compared to that on free protons is the product 
of number abundance times capture rate. Figure |lj 
shows the time evolution of the number abundances for 
free neutrons, protons, a-particles and heavy nuclei, cal- 
culated for t h e sam e stellar trajectory (obtained from 
Liebendorfer , 2002 ) for which the capture rates have 
been evaluated under the assumption of nuclear statis- 
tical equilibrium (NSE). (We note that the commonly 



adopted equations of states (Lattimcr and Swesty, 1991 
Shcn et ai, 1998a| ) yield somewhat larger Y p fractions 



than obtained in NSE.) Importantly the number abun- 
dance of heavy nuclei is significantly larger than that of 
free protons (by more than two orders of magnitude at 
the example point discussed above) to compensate for 
the smaller capture rates on heavy nuclei. It appears 
thus that electron capture on nuclei cannot be neglected 
during the collapse. We note that the average energies of 
the neutrinos produced by capture on nuclei are signifi- 
cantly smaller than for capture on free protons making 
this process a potentially important source for low-energy 
neutrinos. 

The neutron number N = 40 is not magic in nuclear 
structure, nor for stellar electron capture rates. Thus the 
anticipated strong reduction of the capture rate on nu- 
clei will not occur and we expect capture on nuclei to be 
an important neutronization process probably until neu- 
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trino trapping. The magic neutron number N = 50 is 
also no barrier as for nuclei like 93 Kr (N = 57), the neu- 
tron j?/-shell is nearly completely occupied, but due to 
correlations protons occupy, for example, the gg/ 2 orbital, 
and thus unblock GT transitions by allowing transforma- 
tions into gg/2, 9t/2 neutrons. The description of electron 
capture on nuclei in the collapse simulations needs to be 
improved. 

In the current simulations the inverse reaction rates of 
the weak processes listed above are derived by detailed 
balance. Thus an improved description of electron cap- 
ture will then also affect the neutrino absorption on nu- 
clei, although this process is strongly suppressed by Pauli 
blocking in the final state. 



B. Neutrino rates 
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FIG. 15 Normalized neutrino spectra for stellar electron cap- 
ture on the six most important 'electron-capturing nuclei' in 
t he presupernova m odel of a 15 Mq star, as identified in 
The stellar parameters are T = 7.2 x 10 9 
" 3 and Y e =0.43. The solid lines rep- 



( |Hcger et al.[ pOOle 
K, p= 9.1 x 10 a g cm 
resent the spectra derived from the shell model electron cap- 
ture rates. The dashed and dashed- dotted lines cor r espond 
to p aramctrization s recommended by Langanke et al. (20011) 
and Bruenn ( 1981: ), respectively. 



In the capture process on nuclei, the electron has to 
overcome the Q-values of the nuclei and the internal ex- 
citation energy of the GT states in the daughter, so the 
final neutrino energies are noticeably smaller than for 
capture on free protons. Typical neutrino spectra for 
a presupernova model are shown in figure [l5|. In this 
stage of the collapse the neutrino energies are sufficiently 
small that they only excite allowed transitions. Conse- 
quently neutrino-nucleus cross sections for p/-shell nuclei 
can be determined on the basis of GT distributions de- 
termined in the shell model. In the later stage of the 
collapse, the increased density results also in higher en- 
ergy electrons (E e ~ p 1 ^ 3 ) which in turn, if captured by 
protons or nuclei, produce neutrinos with energies larger 
tha n 15-20 MeV. For such neutrinos forbidden (mainly 



A = 1 dipole and spin-dipole) transitions can significantly 
contribute to the neutrino-nucleus cross section. Such 
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FIG. 16 Differential 56 Fe(z/ e , e") 56 Co cross section for the 
KARMEN neutrino spectrum, coming from the decay at rest 
of the muon, as a function of the excitation energy in 56 Co. 
The figure shows the allowed contributions, while the inset 
gives the contribution of the 1~ and 2~ multipolarities. 



a situation is shown in figure 16 which shows the dif- 
ferential cross section for the process 56 Fe(^ e , e~) 56 Co 
computed using the Shell-Model for the Gamow- Teller 
contri bution and th e CRPA for the forbidden contribu- 
tions (Kolbe et al., 1999a). The calculation adopts a 



neutrino spectrum corresponding to a muon decaying 
at rest. The average energy, E v w 30 MeV, and mo- 
mentum transfer, q w 50 MeV, represent the maximum 
values for v e neutrinos expected during the supernova 
collapse phase, i.e., the maximum contribution expected 
from forbidden transitions to the total neutrino-nucleus 
cross sections. In this particular case the 1 + multipole 
(Gamow- Teller at the q — limit) represents 50% of 
the cross section. The 56 Fe(^ e , e~) 56 Co cross section for 
neutrinos from muon-decay at rest has been measured 
by the Karmen collaboration. The measured cross sec- 
tion (2.56±1.08(stat)±0.43(syst) x 10~ 40 cm 2 ) JZcitnitzj 



1994) agrees with the result calculated in the shell model 
(allowed transitions) plus CR PA (forbidden trans itions) 
approach (2.38 x 10" 40 cm 2 ) flKolbe et all Il999a|). 



Although the most important neutrino reactions dur- 
ing collapse are coherent elastic scattering on nuclei 
and inelastic scattering off electrons, it has been noted 
(Bruenn and Haxton, 1991; Haxton, 1988Q that neutrino- 
induced reactions on nuclei can happen as well. Us- 
ing 56 Fe as a representative nucleus, Bruenn and Hax- 
ton concluded that charged-current {v e ,e~) reactions do 
not have an appreciable effect on the evolution of the 
core during infall, due to the high-threshold for neutrino 
absorption. Based on shel l model calcula t ions of the 
GT strength distributions, [Bampaio et al. (2001) con- 
firmed this finding for other, more relevant nuclei in 
the core composition. The same authors showed that 
finite-temperature effects can increase the (v e , e~) cross 
sections for low neutrino energies drastically ( ^ampaio 



et al. , 2001 ). But this increase is found to be significantly 



smaller than the reduction of the cross section caused by 
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Pauli blocking of the final phase space, i.e. due to the 
increasing electron chemical potential. This environmen- 
tal effect ensures that neutrino absorption on nuclei is 
unimportant during the collapse compared with inelastic 
neutrino-electron scattering. 
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FIG. 17 Inelastic neutrino cross sections for 56 Fe (left) and 
59 Co (right) as function of initial neutrino energy and for 
selected temperatures (upper part). Only allowed Gamow- 
Teller transitions have been considered. Temperatures are in 
MeV. For T = 0, the cross section is calculated for the ground 
state only. At T > 0, the cross sections have been evaluated 
for a thermal ensemble of initial states. The corresponding 
neutrino energy distribution in the final state is shown in the 
lower part, assuming an initial neutrino energy of E v — 7.5 
MeV. Due to threshold effects a significant portion of the 
neutrinos are upscattered in energy for even-even nuclei. 



Brucnn and Haxton ( 1991 ) observed that inelastic neu- 
trino scattering off nuclei plays the same important role 
of equilibrating electron neutrinos with matter during in- 
fall as neutrino-electron scattering. The influence of fi- 
nite temperature on inelastic neutrino-nucleus scattering 
was studied in (Fuller and Meyer, 1991), us ing an inde- 
pendent parti cle model. While the study in (Brucnn and 



Haxton 



1991 ) was restricted to 56 Fe, additional cross sec- 
tions have been calculated for inelastic scattering of neu- 
trinos on other nuclei based on modern shell-model GT 
strength distributions ( [Sampaio et al. , 2002b). Again, 



for low neutrino energies the cross sections are enhanced 
at finite temperatures (figure 17). This is caused by the 
possibility that, at finite temperatures, the initial nucleus 
can reside in excited states which can be connected with 
the ground state by sizable GT matrix elements. These 
states can then be deexcited in inelastic neutrino scat- 
tering. Note that in this case the final neutrino energy is 
larger than the initial (see figure [l?]) so that the deexci- 
tation occurs additionally with larger phase space. Un- 
til neutrino trapping there is little phase space blocking 
i n ine lastic neutrino-nucleus scattering. Toivancn et al. 
( [2001 ) presented the charged- and neutral-current cross 
sections for neutrino-induced reactions on the iron iso- 



topes 



52-60 



Fe, using a combination of shell model and 



RPA approach. Other possible neutrino processes, e.g. 
nuclear deexcitation by neutrino pair production (BH), 



estimated rates are probably too small for these processes 
to be important during the collapse. 

Finally we remark that coherent elastic scattering on 
nuclei scales like ~ E% so that neutrinos with low energies 
are the last to be trapped. In order to fill this important 
sink for entropy and energy, processes which affect the 
production of neutrinos with low energies can be quite 
relevant for the collapse. Inelastic neutrino scattering on 
nuclei, i ncluding finite tempe r ature effects, is one such 
process (Bruenn and Haxton, 1991). The significantly 
lower energies of the neutrinos generated by electron cap- 
ture on nuclei than the ones generated by capture on free 
protons is another reason to implement these processes 
with appropriate care in collapse simulations. 



C. Delayed supernova mechanism 
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FIG. 18 Sketch of the stellar core during the shock revival 
phase. R v is the neutrinosphere radius, from which neutrinos 
are expected to stream out freely, R ns is the radius of the pro- 
toneutron star, R g the gain radius (see text) and R a the radius 
at which the shock is stalled. The shock expansion is impeded 
by mass infall at a rate M, but supported by convective en- 
ergy transport from the region of strongest neutrino heating 
to the stalled shock. Convection inside the protoneutron star 
(PNS) as well as correlations in the dense nuclear medium 



incre ase the neutrino luminosity (adapted from Janka et al. 
200l|). 



In the delayed supernova mechanism the fate of the 
explosion is determined by several distinct neutrino pro- 
cesses. When the shock reaches the v e neutrinosphere, 
from which v e are expected to stream out freely, elec- 
tron capture on the shock-heated and shock-dissociated 
matter increases the v e production rate significantly. Ad- 
ditionally neutrinos are produced by the transforma- 
tion of electron-positron pairs into vv pairs (equation 
Bjl) . This process is s trongl y temperature dependent (e.g. 



Soyeur and Brown, 1979] ) and occurs most effectively 



in the shock-heated regions of the proto-neutron star. 
Electron-positron pair annihilation and nucleon-nucleon 
bremsstrahlung (equation Q) generate pairs of all three 
neutrino flavors with the same probability and thus are 
the main mechanisms for the production of Vn . v T neutri- 



have been discussed in (Fuller and Meyer, 1991), but the nos and antineutrinos (Hanhart et at, 2001; Hannestad 
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2001; lHanncstad and Raffelt, 1998; Raffclt, 2001 


Sto- 


ica tnd Horvath|, |2002 


: [Thompson and Burrows , 


200 ll 


Thompson et al, 200C 


). The emitted v e and v e neutri- 



nos, however, can be absorbed again by the free nucleons 
behind the shock. Due to the temperature and density 
dependences of the neutrino processes involved, neutrino 
emission wins over neutrino absorption in a region inside 
a certain radius (the gain radius), while outside the gain 
radius matter is heated by neutrino interactions that are 
dominated by absorption of electron neutrinos and an- 
tineutrinos on free nucleons which have been previously 
liberated by dissociation due to the shock (see figure 18). 
As a net effect, neutrinos transport energy across the gain 
radius to the layers behind the shock. Due to the smaller 
abundances, neutrino-induced reactions on finite nuclei 
are expected to contribute only modestly to the shock 
revival. It has been also sugg ested that the shock revival 
is supported by 'preheating' (Haxton, 1988). In this sce- 
nario the electron neutrinos, which have been trapped 
during the final collapse phase and are liberated in a very 
short burst (with luminosities of a few 10 53 erg s _1 lasting 
for about 10 ms), can partly dissociate the matter (e.g. 
iron and silicon isotopes) prior to the shock arrival. As re- 
liable neutrino-induced cross sections on nuclei have not 
been available until recently, the neutrino-nucleus reac- 
tions have not been included in collapse and post-bounce 
simulations. 

To describe the important neutrino-nucleon processes, 
most core collapse simulations use the same lowest or- 
der cros s section for both neu trinos and antineutrinos 
(Brucnn, 1985) ; |Horowitz , 2002), i.e., they neglect terms 
of order E v /M, where E v is the neutrino energy and 
M the nucleon mass. The most important corrections 
to the cross section at this order are the nucleon recoil 
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FIG. 19 Shock trajectories of a 20 Mq star, calculated with 
(dashed) and without (solid) an isoscalar strange axialform 
factor in the neutrino-nucleon elastic cross sections (courtesy 
of M. Liebendorfer). 



and assuming axial-vector dominance, i.e. the cross sec- 
tion scales like a ~ \gA T ~ 9a\ 2 ' a non- vanishing strange 
axial- vector form factor would reduce the elastic scatter- 
ing cross sect ion on neutrons and increase the up elasti c 
cross section (Garvey et al., 1993a, 1992; Horowitz, 2002). 
As the matter behind the shock is neutron-rich, the net 
effect will be a reduction of the neutrino-nucleon elas- 
tic cross section. This increases the energy transfer to 
the stalled shock, however, a simulation has shown that 
this increase is not strong enough for a successful shock 
revival (Liebendorfer et al, 2002, see figure [Tgj) . 

The physics involved in the attempt to revive the shock 
by neutrino heating i s exhaustively reviewed b y lanka 



in Eql (|4a|) ( [Horowitz] , ^UL)2|T ! the recoil correction is 
the same for neutrinos and antineutrinos and decreases 
the cross sections. However, the weak magnetism cor- 
rects the cross sections via its parity-violating interfer- 
ence "witliTniFlIomln^^ 

interference is constructive tor neutrinos and destructive 
for antineutrinos, inclusion of the weak magnetism cor- 
rection increases the neutrino cross section, while it de- 
creases the P-nucleon cross sections. It is then expected 
that corrections up to order E v jM decrease the antineu- 
trino cross section noticeably (by about 25% for 40 MeV 
antineutrinos), while the i/- nucleon cross sections are only 
affect ed by a few percents for E v < 100 MeV (Horowitz, 
200^ ). 

Neutral-current processes are sensitive to possible 
strange quark contributions in the nucleon which would 
give rise to an isoscalar piece g s A in the axial- vector form 
factor besides the standard isovector form factor gA~r 
( [Beise and McKeown| , |1991| ; |Jaffe and Manohaij , [L990[ ), 
The current knowledge on g s A comes from a vp clastic 



andlthc^cak magnetism related to the form factor E 2 e t al\ (§00f) (see also [Burrows and Goshyj [19931). These 



authors show that the fate of the stalled shock does not 
only depend on the neutrino heating above the gain ra- 
dius, but is also influenced by the energy loss in the cool- 



ing region below t h e gain radius (see also Bethc and Wil , 
1985| ). lanka ( 2001 ) also demonstrates the existence 



son 



of a critical value for the neutrino luminosity from the 
neutron star needed to revive the shock. This critical 
luminosity depends on the neutron star mass and radius 
and on the mass infall to the shock. One expects that the 
shock expansion is eased for high mass infall rates, which 
increase the matter pile-up on the neutron star and push 
the shock outwards, and for high v e and v e luminosities 
from the neutron star, which lead to an enhancement of 
neutrino absorption relative to neutrino emission in the 
gain region. 

The explosion depends crucially on the effectiveness 
by which energy is transported by neutrinos to the re- 
gion where the shock has stalled. As stressed before, 
one-dimensional models including so p histic a ted neutrino 



transp ort (e .g. Liebendorfer et al., [2001 ; Rampp and 



i lie cmicm kiiowiuu^c un y A comes uom a up ciasnc 

scat tering experiment performed at Brookhaven yielding Janka , 2000D fail to explode. However, the neutrino en- 
g% = -0.15±0.08 ([Ahrens et al\ , [19871) , but is considered 
rather uncertain (Garvey et al, 1993b| ). With gA = 1-26 



ergy transport is very sensitive to: i) the effect of nucleon- 
nucleon correlations on the neutrino opacities in dense 
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matter and ii) convection both in t he neutrino-heate d 
region and in the proto- neutron star (Janka et al, [2001 ). 

In his pioneering work, [Sawyer (1989) calculated the 
neutrino mean free path, or equivalently the neutrino 
opacity, in uniform nuclear matter and showed that ef- 
fects due to strong interaction between nucleons are im 



por tant. The sam e conclusion has been reached by |Raf- 
felt ~ 



et al. 



( 1996 ) who demonstrated that the average 
neutrino-nucleon cross section in the medium is reduced 
due to spin fluctuations induced by the spin-dependent 
interaction among nucleons. (For earlier calculations of 
the neutrino mean free path in uniform nuclear matter, 



see Friman and Maxwell [1979|; |Iwamoto and Pethick 



1982) ; pawyeij |1975| ). |5awyer| ( |1989| ) exploited the rela- 
tion between the equation of state (EOS) of the matter 
and the long-wavelength excitations of the system to cal- 
culate the weak interaction rates. However, consistency 
between the EOS and the neutrino opacities are more dif- 
ficult to achieve for large energy (qo) and momentum (q) 
transfer of the neutrinos. Here, particle- hole and particle- 
particle interactions are examples of effects which might 
influence the EOS and the neutrino opacities. For the 
following discussion it is quite illuminating to realize the 
similarity of the neutrino-induced excitations of nuclear 
matter with the physics of multipole giant resonances in 
finite nuclei. 

For muon and tau neutrinos, neutral current reactions 
are the only source of opacities. Here, the energy and 
momentum transfer is limited by the matter tempera- 
ture alone. For electron neutrinos the mean free path 
is dominated by charged-current reactions, for which the 
energy transfer is typically of the order of the difference 
between neutron and proton chemical potentials. During 
the early deleptonization epoch of the proto-neutron star 
the typical neutron momenta are large (~ 100-200 MeV) 
and the mismatch of proton, neutron and electron Fermi 
momenta can be overcome by the neutrino momenta. 
This is not longer possible in later stages when the neu- 
trino energies are of order kbT; then momentum conser- 
vation restricts the available phase space for the absorp- 
tion reaction. Pauli blocking of the lepton in the final 
state increases the mean free path for charged-current 
and neutral-current reactions. 

We note an important and quite general consequence 
of the fact that muon and tau neutrinos react with the 
proto-neutron star matter only by neutral-current reac- 
tions: The four neutrino types have similar spectra. Due 
to universality, i/^, v T and v^, v T have identical spec- 
tra. It is usually even assumed that neutrinos and an- 
tineutrinos have the same spectra (one therefore refers 
to the 4 neutrino types unifyingly as v x neutrinos) ex- 
ploiting axial-vector dominance in the neutrino cross sec- 
tions. However, the interference of the axial-vector and 
the weak magnetism components makes the v x spectra 
slightly hotter than the v x spectra. The v x neutrinos 
decouple deepest in the star, i.e., at a higher temper- 
ature, than electron neutrinos and antineutrinos, and 
hence have higher energies. As the matter in the proto- 



neutron star is neutron-rich, electron neutrinos, which 
are absorbed by neutrons, decouple at a larger radius 
than their antiparticles, which interact with protons by 
charged-current reactions. As a consequence decoupled 
electron neutrinos have, on average, smaller energies than 
electron antineutrinos. Calculated supernova neutrino 



spectra can be found in (Janka and Hillebrandt, 1989 
Yamada et al. , 1999| ), which yield the average energies 



of the various supernova neutrinos approximately as 
(E Ue ) = 11 MeV, {E 9 „) = 16 MeV, and (E Vx ) = 25 
MeV. |Burrows et ~al\ ( p000| ) find the same hierarchy, but 
somewhat smaller average neutrino energies. 

For a much deeper and detailed description of the neu- 
tri no mean free paths in dense matter t he rea der is refered 
to ( Prakash et al , 2001 ; Reddy et al[ 1999 ) and the ear- 



lier work (Burrows and Sawyer, 1998, 1999] ; Reddy et al 
1998 ). We will here only briefly summarize the essence 



of the work presented in these references. 
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FIG. 20 Ratio of neutrino mean free paths in neutron matter 
calculated in RPA and Hartree-Fock approaches at various 
temperatures ( Margueron et al. , 2002 ) . The interaction is the 
Gogny force DIP. The neutrino energy is taken as E v — 3T. 



Collapse simulations describe neutrino opacities typi- 
cally on the mean-field level or even by a nucleon gas. 
Then an analytical expression can be derived for the vec- 
tor and axial-vector response of the medium which in 
turn determines the charged- and neutral-current cross 
sections. Effects due to the strong interaction between 
nucleons are considered by a medium-dependent effec- 
tive mass in the dispersion relation. Like in finite nu- 
clei, collective excitations in nuclear matter arise due to 
nucleon-nucleon correlations beyond the mean-field ap- 
proximation. As it is believed that single-pair excitations 
dominate over multi-pair excitations for the kinematics 
of interest to neutrino scattering and absorption, it ap- 
pears to be sufficient to determine the vector and axial- 
vector response, in a first step, within the Random Phase 
Approximation (RPA). Assuming that the interaction is 
short-ranged compared to the wavelength of the excita- 
tions, it is justified to retain only s-wave components in 
the interaction which in turn can be related to Fermi- 
liquid parameters. It is found that the repulsive nature of 
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the parameter G' , which is related to the isovector spin- 
flip or giant Gamow- Teller resonances in nuclei, induces 
a collective state in the region u>/q ~ Vp (vf is the Fermi 
velocity) , while the cross section is reduced at smaller en- 
ergies. However, these smaller energies are important for 
the neutrino mean free path at nuclear matter densities 
(po) or smaller densities. Assuming a typical neutrino 
energy E v « 3T (corresponding to a Fermi-Dirac distri- 
bution with temperature T and zero chemical potential) 
RPA correlations increase the neutral-current neutrino 
mean-free path (see figure p0[ ) at low temperatures and 
for p = po, compared with the mean- field result. An en- 
hancement due to RPA correlations is also found for neu- 
trino absorption mean-free paths for neutrino-trapped 



1998). Such a damping is possible in regions where neu- 
trinos are still strongly coupled to matter, but neutrino 
opacities are not too high t o make neutrino transport 

~ ( |1996| ) the con- 



insufficient. In the model of Keil et 



mat ;er. Like in the case of neutrino-induced reactions on Dineva, 1996) 



vective mixing occurs very deep inside the core where the 
neutrino opacities are high; no damping of the convection 
by neutrinos is then found. 

Wilson and Mayle attempted to simulate convection in 
their spherical model by intro ducing neutron-fing e rs an d 
found successful explosions ( Wilson and Mayle , 1993 ). 
This idea is based on the assumption that energy trans- 
port (by three neutrino flavors) is more efficient than 
lepton number transport (only by electro n neutrinos). 
However, this assumption is under debate (Bruenn and 



finite nuclei (see above), finite-temperature effects allow 
that nuclear excitation energy is transferred to the neu- 
trino in inelastic scattering processes. This contributes 
to the cooling of the nuclear matter and increases the 
neutrino energy in the final state. 

Neutrino heating is maximal in the layer just above the 
gain radius. The energy transport from this region to the 
shock, which is stalled further out, can be supported by 
convective overturn and might lead to successful explo- 
sions, as has been demonstrated in several simulations 
with two-dimensional hydrodynamics treatme nt of the 
region between the gain radius and the shock ( Burrows 
et afj, |1995| ; |Hcrant et al\ |1994| ; [Janka and Mulleij , fl996j) 
The effect of convection is twofold (Janka et al[ pOOl 



At first, heated matter is transported outwards to cooler 
regions where the energy loss due to neutrino emission 
is reduced (the neutrino production rate for electron and 
positron captures on nucleons depends strongly on tem- 
perature). Second, cooler matter is brought down closer 
to the gain radius where the neutrino fluxes are larger and 
hence the heating is more effective. While this picture is 
certainly appealing, it is not yet clear whether multi- 
dimensional simulations will indeed lead to explosions as 
the two-dimensional studies did not include state-of-the- 
art Boltzmann neutrino transport, but treated neutrino 
transport in an approximate manner. In fact, in a sim- 
ulation with an improved treatment of neutrino trans- 
port the convective overturn was fou nd not to be strong 



enou gh to revive the stalled shock (Mezzacappa et al 



1998 ) 



The shock revival can also be supported by convection 
occuring inside the protoneutron star where it is mainly 
driven by the negative lepton gradient which is estab- 
lished by the rapid loss of leptons in the reg i on ar ound 



the neutrinosphcrc (Burrows, 1987; Keil et al, 1996). By 



this mode, lepton-rich matter will be transported from 
inside the protoneutron star to the neutrinosphere which 
increases the neutrino luminosity and thus is expected to 
help the explosion. The simulation of protoneutron star 
convection is complicated by the fact that neutrinos and 
matter are strongly coupled. In fact, two-dimensional 
simulations found that neutrino transport can equilibrate 



VI. NUCLEOSYNTHESIS BEYOND IRON 

While the elements lighter than mass number A ~ 60 
are made by charged-particle fusion reactions, the heav- 
ier nuclei are made by neutron captures, which have to 



compete with 8 decays. Already (Burbidge et al, 1957 



otherwise convective fluid elements (Mezzacappa et al 



Cameron, 1957) realized that two distinct processes are 
required to make the heavier elements. This is the slow 
neutron-capture process (s-process), for which the 8 life- 
times Tp are shorter than the competing neutron capture 
times t„. This requirement ensures that the s-process 
runs through nuclei in the valley of stability. The rapid 
neutron-capture process (r-process) requires r n -C 7/3 . 
This is achieved in extremely neutron-rich environment, 
as r n is inversely proportional to the neutron density 
of the environment. The r-process runs through very 
neutron-rich, unstable nuclei, which are far-off stability 
and whose physical properties are often experimentally 
unknown. 

Weak-interaction processes play interesting, but dif- 
ferent roles in these processes. The half-lives of the 
/3-unstable nuclei along the s-process path are usually 
known with good precision. However, the nuclear half-life 
in the stellar environment can change due to the thermal 
population of excited states in the parent nucleus. This 
is particularly interesting if the effective lifetime is then 
comparable to t„, leading to branchings in the s-process 
path, from which the temperature and neutron density 
of the environment can be determined. Several recent 
examples are discussed in the next subsection. For the r- 
process, /3-decays are probably even more crucial. They 
regulate the flow to larger charge numbers and determine 
the resulting abundance pattern and duration of the pro- 
cess. Except for a few key nuclei, /3-decays of r-process 
nuclei have to be modelled theoretically; we will briefly 
summarize the recent progress below. Although the r- 
process site is not yet fully determined, it is conceivable 
that it occurs in the presence of extreme neutrino fluxes. 
As we will discuss, neutrino-nucleus reactions can have 
interesting effects during and after the r-process, perhaps 
allowing for clues to ultimately identify the r-process site. 
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A. S-process 



The analysis of the solar abundances have indicated 
that two components of the s-process have contributed 
to the synthesis of elements heavier than A ~ 60. The 
weak component produces the elements with A < 90. Its 
site is related with helium c ore burning of CNO mate- 



rial in more massive stars (Couch et al, 1974; Kappeler 



et a (j, |199^ ). The main component, which is responsible 



for the heavier s-process nuclides up to Pb and Bi, is asso- 
ciated with helium flashes occurring during s hell burning 



in low-mass (asymptotic giant branch) stars (Busso et 
1999| ). The 22 Ne(a,n) 25 Mg and 13 C(a,n)">0 reactions 



are believed to be the supplier of neutrons for the weak 
and main components, respectively. The s-process abun- 
dances N s are found to be inversely proportional to the 
respective (temperature averaged) neutron capture cross 
sections 



as ex pected for a steady-flow picture (Bur 
bidg e et a/.| , 1957j ), which, however, breaks down for the 



cxtcmcly small cross sections at the magic neutron num- 
bers. As a consequence, the product N s ■ (a) exhibits 
almost constant plateaus between the magi c neutron 
numbers, separated by pronounced steps, (e.g. Kappelei , 



1999) 
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FIG. 21 The s-process reaction path in the Nd-Pm-Sm re- 
gion with the branchings at A = 147, 148, and 149. Note 
that 148 Sm an d 150 Sm ar e shiel ded against r-process /3-decays 
(adapted from Kappelei, 1999). 



The neutron density of the stellar environment during 
the main s-process component can be determined from 
branching points occuring, for example, in the A = 147- 
149 mass region (see figure ^l]). Here the relative abun- 
dances of the two s-process-only isotopes 148 Sm and 
i50g m ^2 = 62), which are shielded against r-process 
contributions by the two stable Nd (Z = 60) isotopes 
148 Nd and 150 Nd, is strongly affected by branchings oc- 
curing at 147 Nd and, more importantly, at 148 Pm and 
147 Pm. As the neutron captures on these branching 



at 



nuclei will bypass 148 Sm in the flow pattern, the 150 Sm 
N s (cr) value will be larger for this nucleus than for 148 Sm. 
Furthermore, the neutron capture rate A„ is proportional 
to the neutron density N n . Thus, N n can be determined 
from the relative 150 Sm/ 148 Sm abundances, r e sultin g in 
N n = (4.1 ± 0.6) x 10 8 cm" 3 ( |Toukan et aJ\ , |1995| ). A 
similar analysis for the weak s-process component yields 



neutron densities of order (0.5-1.3) x 10 s cm 3 (Wal- 
1986a ,jj|). We stress that a 10% determina- 



ter et al. 



tion of the neutron density requires the knowledge of 
the involved neutron capture cross sections with about 
1% accuracy (Kappeler et al., 1998), which has yet not 



been achieved for unstable nuclei. Improvements are ex- 
pected from new time-of-flight facilities like LANSCE at 
Los Alamos or NTOF at CERN. 
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region (solid lines). For 176 Lu the ground state and isomer 
are shown separately. Note that 176 Lu an d 176 Hf a re shielded 



against r-process /3-decays (adapted from Doll et al., 1999) 
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FIG. 23 The figure shows an schematic energy-level diagram 
of 180 Ta and its daughters illustrating the possibility for ther- 
mally enhanced decay of 180 Ta m in the stellar environment of 
the s-process. The inset sho ws the photon den sity at s-process 
temperature (adapted from Belie et al., 1999). 



25 



The temperature of the s-process environment can be 
'measured', if the j3 half- life of a branching point nu- 
cleus is very sensitive to the thermal population of ex- 
citcd n uclea r levels. A prominent example is 176 Lu (Beer 
et a|I |1981 ). For mass number A — 176, the /3-decays 
from the r-process terminate at 176 Yb (Z = 70), mak- 
ing 176 Hf (Z = 72) and 176 Lu (Z = 71) s-only nu- 
clides (see figure |2^) . Besides the long-lived ground state 
(t 1/2 = 4.00(22) x 10 10 y), 176 Lu has an isomeric state 
at an excitation energy of 123 keV (tx/2 = 3.664(19) h). 
Both states can be populated by 175 Lu(n, 7) with known 
partial cross sections. At 176 Lu, the s-process matter 
flow is determined by the competition of neutron cap- 
ture on the ground state and /3-decay of the isomer. 
But importantly, the ground and isomeric states cou- 
ple in the stellar photon bath via the excitation of an 
intermediate state at 838 keV (figure ^3] shows a simi- 
lar process taking place in 180 Ta), leading to a matter 
flow from the isomer to the ground state, which is very 
temperature-dependent. A recent analysis of the 176 Lu 
s-process branching yields an environment temperature 
of T = (2.5-3.5) x 10 s K (poll et al\, ggggh. 



Similar finite-temperature effects play also an impor- 
tant role in the s-process production of 180 Ta. This is 
the rarest isotope (0.012%) of nature's rarest element. 
It only exists in a long-lived isomer (J 1 * = 9~) at an 
excitation energy of 75.3 keV and with a half-life of 
ti/2 > 1-2 x 10 15 y. The 1+ ground state decays with 
a half-life of 8.152(6) h, mainly by electron capture to 
180 Hf. While potential s-process production paths of 
the 180 Ta isomer have been pointed out, the survival of 
this state in a finite-temperature environment has long 
been questionable. While a direct electromagnetic de- 
cay to the ground state is strongly suppressed due to 
angular momentum mismatch, the isomer can decay via 
thermal population of intermediate states with branch- 
ings to the ground state (see figure |23| ). By measur- 
ing the relevant electromagnetic coupling strength, the 
temperature-dependent half-life, and thus the 180 Ta sur- 



vival rate, has been determined by Belie et al. (199E) 



under s-process conditions as t\i% < 1 y, i.e. more than 
15 orders of magnitude smaller than the half-life of the 
isomer! Accompanied by progress in stellar modeling of 
the convective modes during the main s-process com- 
ponent (which brings freshly produced 180 Ta to cooler 
zones, where it can survive more easily, on timescales of 
days) it appears now likely that 180 Ta is partly mad e 



within s-process nucleosynthesis (Wisshak et al. , |2001 ) 



The p- proc ess (Rayet et al. , 1995 ) and neutrino nucle- 
osynthesis ( Wooslcy et al. , 1990|) have been proposed as 
alternative sites for the 180 Ta production. 

The two neighboring isotopes 186 Os and 187 Os are s- 
process only nuclides, shielded against the r-process by 
186 W and 187 Rh. The nucleus 187 Rh has a half-life of 42 
Gy which is comparable with the age of the universe. As 
the 187 Rh decay has contributed to the observed 187 0s 
abundance, the Os/Rh abundance ratio can serve as a 
sensitive clock for the age of the universe, once the s- 



process component is subtracted from the 187 0s abun- 
dance ( Clayton , 1964 ). To determine the latter precise 
measurements of the neutron capture cross sections on 
186 0s and 187 0s are required, which are in progress at 
CERN's NTOF facility. A potential complication arises 
from the fact that 187 0s has a low-lying state at 9.8 
keV which is populated equally with the ground state 
at s-process temperatures. The respective neutron cap- 
ture cross sections on the excited state can be indi- 
rectly determined from (n,n') measurements. Further- 
more one has to consider that the half-life of 187 Rh is 
strongly temperature-dependent: If stripped of its elec- 
trons, the half-life is reduced by 9 orders of magnitude 
to t-i/ 2 — 32.9 ± 2 y as measured at the GSI storage 



ring ( Bosch et al. , 1996 ). However, the GSI data can 
be translated into a log ft value from which the 187 Rh 
half-life can be deduced for every ionization state. 

The decay of totally ionized 187 Rh is an example for 
a bound state /3-decay, i.e., the decay of bare 187 Rh 75+ 
to continuum states of 187 0s 76+ is energetically forbid- 
den, but it is possible if the decay electron is captured in 
t he K-shcll (Q fl =73 keV ) or in the L-shell ((3/3=9.1 keV) 
(Johnson and Sofl, 1985). We note that the decay of neu- 
tral i8Y Rh is energetically allowed for a first-forbidden 
transition to the 187 Os ground state; the long half-life re- 
sults from the small matrix element and Qp value related 
to this transition. Another example of a bound state /3- 
decay with importance for the s-process occurs for 163 Dy. 
This nucleus is stable as a neutral atom, but, if fully ion- 
ized, can decay to 163 Ho with a half-life of 47^ d; the 
measurement of this half-life at the GSI stor age ring was 



the first observation of bound state /3-decay ( Jung et al 



1992). The consideration of the 163 Dy decay at s-process 



conditions has been found to be essential to explain the 
abundance of the s-process only nuclide 164 Er, which is 
produced by neutron capture on 163 Ho, the daughter of 
the bound state /3-decay of 163 Dy, and the subsequent 
/3-decay of 164 Ho to 164 Er. 



B. R-process 

Phcnomcnological parameter studies indicate that the 
r-process occurs at temperatures around T ~ 100 keV 
and at extreme neutron fl uxes (neutron numb er densi- 
ties n > 10 20 cm~ 3 ) (e.g. |Cowan et al\, |l99l[). It has 



also been demonstrated that not all r-process nuclides 
can be made simultaneously at the same astrophysical 
conditions (constant temperature, neutron density), i.e. 
the r-process is a dynamical p rocess with changin g (dif- 



ferent) conditions and paths ( Kratz et al , 1993| ). In a 



good approximation, the neutron captures proceed in 
(n, 7) (7, n) equilibrium, fixing the reaction paths 
at rat her low neutron sep aration energies of S n ~ 2-3 
MeV ( |Cowan et al. , 1991 ) , implying paths through very 
neutron-rich nuclei in the nuclear chart as is shown in 
figure ^J. While the general picture of the r-process ap- 
pears to be well accepted, its astrophysical site is still 
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open. The extreme neutron fluxes point to explosive sce- 
narios and, in fact, the neutrino-driven wind above the 
nascent neutron star in a core-collapse super no va is the 



curren tly fa vored model ( Witti et al , 1994t]]b| |Wooslcy 



et ofl 1994). But also shock-pro cessed helium shells in 



type II s upcrnovac ( Iruran et al , [2001 ) and neutron sta r 
mergers ( Freiburghaus et al, 1999b; Rosswog et al, 200C) 
are investigated as possible r-process sites. Recent mc- 
teoritic clues (Qian et al, 1998] ; Wasscrburg et al, 1996) 
as well as observations of r-proccss abundances in low- 
metallicity stars ( (Bneden et al , 200C ) point to more than 



one distinct site for the solar r-process nuclides. 



In a n important astronomical observation Cayrel et al. 
( [2001 ) have recently detected the r-process nuclides tho- 



rium and uranium in an old galactical halo star. As these 
two nuclei have halflives, which are comparable to the ex- 
pected age of the universe, their measured abundance 
ratio serves as a sensitive clock to determine a lower 
limit to this age, provided their initial r-process produc- 
tion abundance ratio can be calculated with sufficient 



As relevant nuclear input, r-process simulations require 
neutron separation energies (i.e. masses), half- lives and 
neutron capture cross sections of the very neutron-rich 
nuclei on the various dynamical r-process paths. As cur- 
rently only few experimental data are known for r-process 
nuclei, these quantities have to be modelled. Tradition- 
ally this has been done by global models which, by fit- 
ting a certain set of parameters to known experimental 
data, are then being used to predict the properties of all 
nuclei in the nuclear landscape. Arguably the most im- 
portant input to r-process simulations are neutron sepa- 
ration energies as they determine, for given temperature 
and neutron density of the astrophysical environment, 
the r-process paths. The most commonly used mass tab- 
ulations are based on the microscopic-macroscopic finite- 



range droplet model (FRDM) approach (Moller et al 



1997 ) or the Extended Thomas-F ermi model with Struti - 
nski Integral (ETFSI) approach (Aboussir et al, 1995) 



accuracy. Cayrel et al (2001) have used their observed 
238 U/ 232 Th abundance ratio from the star CS31082-001 
and the r-process model predictions from (Cowan et al 



199? ]; poriely and Clerbaux , 199E) to deduce a value of 
12.5 ± 3.3 Gyr for the age of the star. Recently, a refined 
analysis of the CS31082-001 spectra has led to a signifi- 
cant improvement in the derived abundances which pro- 



In more recent developments mass tabulations have been 
developed ado pting paramctriz a tions inspired by shell 
model results (Duflo and Zukei, 1995] ) or calculated on 
the basis of nuclear many-body theories like the Hartree- 
Fock model with a BCS treatment of pairing (Goriely 
et al, 2001). Special attention has been paid recently 



also to the 'shell quenching', i.e. the observations made 
in HFB calculations that the shell gap at magic neu- 
tron numbers is less pronounced in v ery neutron-rich nu- 



vides now an age estimate of 14.0 ± 2.4 Gyr (Hill et al, clei t han in nuclei close to stability ( [Dobaczewski et 



2002 ) . The effect of different nuclear physics input on the 
r-process production of U a nd Th h ave been stu d ied b y 
Goriely and Arnould| fl200l| ) and by |Schatz et al\ ( |2002| ) 



1994). Such a vanishing of the shell gap has been exper- 
imentally verified for the magic neutron number N — 20 



(iGuillcmaud- Mueller et al , 1984; Motobayashi et al 



1995). The confirmation of the predicted quenching at 
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the N — 82 shell closure is t he aim of consider able cur- 
rent experimental activities (Kratz et al, [2000 ). Recent 



theoretical s t udies on th is topic are reported in (3harma 
and Farhan, 2001, 2002). The potential importance of 



shell quenching for the r-process rests on the observa- 
tion flChen et al\ , |1995| ; |Pfeiffer et~al\ |1997| ) that it can 
correct the strong trough just below the r-process peaks 
in the calculated r-process abundances encountered with 
conventional mass models. Neutron capture cross sec- 
tions become important, if the r-process flow drops out 
of (n,7) (li n ) equilibrium, which happens close to 
freeze-out when the neutron source ceases. They can also 
be relevant for nuclides with small abundances for which 



no flow equilibrium is built up (Burman and Engel, 2001 



FIG. 24 The r-process occurs under dynamically changing 
astrophysically conditions which affect the reaction pathway. 
The figure shows the range of r-process paths, defined by 
their waiting point nuclei. After decay to stability the abun- 
dance of the r-process progenitors produce the observed solar 
r-process abundance distribution. The r-process paths run 
generally through neutronrich nuclei with experimentally un- 
known masses and halflives. In this calculation a mass formula 
based on the Extended Thomas Fermi model with Strutinski 
Integral (ETFSI) and special treatment of shell quenching 
(see text) has been adopted, (courtesy of Karl-Ludwig Kratz 
and jHeiidrik Sclialz). 



3urman et al, 1997). In the following we will summarize 
recent progress in calculating half-lives for nuclei on the 
r-process paths. 



1. Ha If- Lives 

The nuclear half-lives determine the relative r-process 
abundances. In a simple /3-flow equilibrium picture the 
elemental abundance is proportional to the half-lif e, with 
some corrections for /3-delayed neutron emission (Kratz 
et al, 1988). As r-process half- lives are longest for the 
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magic nuclei, these waiting point nuclei determine the 
minimal r-proccss duration time; i.e. the time needed to 
build up the r-process peak around A ~ 200 via matter- 
flow from the seed nucleus. We note, however, that this 
time depends also crucially on the r-process path. 

Pioneering experiments to measure half-lives of 
neutron-rich isotopes near the r-process path succeeded 
in determining the half- lives of two N — 50 ( 80 Zn, 79 Cu) 
and two A = 82 ( 129 Ag, 130 Cd) waiting point nucle i 



(pill et al 



1986 



Kratz et 



1986 ; Pfeiffer et 



Pfeiffer et al. (2001) reviewed the experimental informa- 
tion on r-process nuclei. These data play crucial roles 
in constraining and testing nuclear models, which are 
still necessary to predict the bulk of half-lives required 
in r-process simulations. It is generally assumed that the 
half-lives are determined by allowed Gamow- Teller (GT) 
transitions. The calculations of /3-decays require usu- 
ally two ingredients: the GT strength distribution in the 
daughter nucleus and the relative energy scale between 
parent and daughter, i.e. their mass difference. However, 
the j3 decays only probe the weak low-energy tail of the 
GT distributio ns. Only a few pe rcent of the 3(N — Z) 



Ikeda sum rule (Ikeda et al., 1963) lie within the Qp win- 



dow (i.e. at energies accessible in /3-decay), the rest being 
located in the region of the Gamow- Teller resonance at 
higher excitation energies. Due to the strong E 5 energy 
dependence of the phase space /3-decay rates are very 
sensitive to the correct description of the detailed low- 
energy GT distribution and its relative energy scale to 
the parent nucleus. This explains why different calcu- 
lations of the /3-decay half-lives present large deviations 
among them. 

Because of the huge number of nuclei relevant for the 
r-process, the estimates of the half-lives are so far based 
on a combination of global mass models and the quasi- 
particle random phase approximation (QRPA), the lat- 
ter to calculate the GT matrix eleme nts. Examples o f 



these models are the FR DM/QRPA flMollcr et al\, |1997| ) 
and the ETFSI/QRPA ( [Borzov and Goriely| , |2000| ). Re- 
cently, calculations based on the self-consistent Hartree- 
Fock-Bogoliubov plus QRPA model became available for 
r-process waiting point nuclei with magic neutron num- 
bers N = 50, 82 and 126 jEngel et cd\ , |L999[ ). The 
presence of a closed neutron shell has allowed also for 
t he study of the s e nu c lei by shell-model calculations 



(Martmez-Pincdo, 2001; Martmcz-Pincdo and Langanke 



1999), which is the method of choice to determine the 
/3-decay matrix elements. However an adequate calcu- 
lation of the nuclear masses f or heavy nuclei is yet pro- 
hibitive in the shell model. In Martmez-Pinedo and Lan- 



were adopted from the global model of Duflo and Zukci 
( |1995|). Figure ^ compares the half-lives predicted by 



ganQ ( |l999 ), and Martmcz-Pinedo| (|2001) the masses 



ciencies. For example, the FRDM/QRPA half- lives show 
a significant odd-even staggering which is not present 
in the data, while the ETFSI/QRPA half- lives appear 
globally too long. The HFB/QRPA and SM approaches 
obtain half-lives in reasonable agreement with the data 
and predict shorter half-lives for the unmeasured waiting 
point nuclei with N = 82 than the global FRDM/QRPA 
and ETFSI/QRPA approaches. For N = 126 there is no 
experimental information so the different models remain 
untested. While HFB calculations for the half-lives of all 
the r-process nuclei are conceivable, similar calculations 
within the shell-model approach are still not feasible due 
to computer memory limitations. 

While the Qp value for the decay of neutron-rich r- 
process nuclei is large, the neutron separation energies 
are small. Hence /3 decay can lead to final states above 
neutron threshold and is accompanied by neutron emis- 
sion. If the r-process proceeds by an (n, 7) <^ (7, n) 
equilibrium the /3-delayed neutron emission probabilities, 
Pp,m only play a role at the end of the r-process when the 
neutron source has ceased and the produced nuclei decay 
to stability. The calculated Pp, n values are very sensitive 
to both the low-energy GT distribution and the neutron 
threshold energies. No model describes currently both 
quantities simultaneously with sufficient accuracy. Fig- 
ure]^ compares the Pp <n computed in the FRDM/QRPA 
and SM approaches for N — 82 and 126. The predictions 
between different models can be quite different in ex- 
perimentally non-determined mass regions. For the SM 
approach the error bars indicate the sensitivity of the 
computed Pp, n values to a change of ±0.5 MeV in the 
neutron separation energies of the daughter nucleus; the 
effect can be large. 

It has been pointed out that first-forbidden transitions 
might have import ant contribution s in some nuclei close 
to magic numbers ( Blomqvist et al. , 1983| ; Homma et al. 



1996| ; [Korgul et alL 2001). A systematic inclusion of first- 
forbidden transitions in the calculation of r-process beta- 
decay half-lives in any of the many-body methods used to 
describe the GT contributions has not been done. How- 



ever, a first attempt towards this goal (Mollcr et al. 



2002) 



has combined first-forbidden transitions estimated in the 
Gross Theory flTakahashi et al. , 1973| ) with GT results 
taken from QRPA calculations. No significant changes 
compared to r-process studies, which consider only the 
GT contributions to the halflives, have been observed 



(|Kratz| , gOOgj ). 

The presence of low-lying isomeric states in r-process 
nuclei might change the effective half-lives in the stellar 
environment. Currently no estimates of these effect exists 
except for odd-A nuclei with N = 82. Here shell-model 
calculations predict half-lives for the iso meric states very 
similar to the ground state half- lives ( Martmcz-Pincdc 



131 In, 1311 Cd and ^Ag ar c know n experimentally ( Kratz 
et a{], |1986| ; pfeiffer et al\ , |2001| ). The comparison of the 



the |diffcrcnt approaches. For N = 82, the half- lives of and Langanke , 1999] ). The halflives of the isomeric state 



predictions of the different models with the few experi- 
mental r-process benchmarks reveals some of their insuffi- 



in lzy Ag has also been estimated within the QRPA ap- 
proach and by a second shell model calculation finding 
values about a factor of 2 larger than the 129 Ag ground 



state half-life (Kratz, 2001). This reference also reports 
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FIG. 25 The left panels show the half-lives of the r-process waiting point nuclei with neutron numbers N — 82 and 126, 
obtaine d within different th eoretical models. In the case of N = 82 the half-lives of 131 In, 130 Cd and 129 Ag are experimentally 



known (Pfeiffer et al 



2001). The right panels compare the beta-delayed neutron emission probabilities for the N = 82 (upper) 
waiting-point nuclei as calculated in two different models. For the shell-model (SM) probabilities the 



and N = 126 (lower 

error bars indicate the sensitivity of the calculations to a change of ±0.5 MeV in the neutron separation energies 



about the first attempt to measure the half-life of the 
isomeric state. 

For heavy nuclei (Z > 84) some final states populated 
by /3-decay in the daughter nucleus can also decay by 
fission ( powan et al. , 1991). The relevant beta-delayed 
fission probabilities de pend sensitively on the modelling 
of t he fission bar riers ( Howard and Mollei , 1980| ; Mam- 



dour et al, 2001) 



Borzov and Goriely have studied the influence of the (3 
half- lives on the r-process abundances within two distinct 
scenarios: the canonical r-process picture with a exposure 
of the seed nucleus 56 Fe by a constant neutron density 
and temperature for a fixed duration time (2.4 s) and 
the neutrino-driven wind model. In the canonical model 
the location of the r-process abundance peaks depends on 
the masses, but not on the j3 half-lives, which act only 
as bottlenecks for the matter-flow to more massive nu- 
clei. In the dynamical neutrino-driven wind model the 
half-lives affect the abundance distribution. This comes 
about as at later times in this model the environmen- 
tal conditions shift the r-process to more neutron-rich 



-floW 



nuclei. Long half-llvcH ihen Imply that ihe matter- 
reacites Llic magic neuti'oii numbers later, i.e. for mure 




FIG. 26 Abundances of r-process nuclides calculated in a dy- 
namical r-process model and for different global sets of /3- 
decay half-lives. In the dynamical r-process the matter flow 
timescale competes with the nuclear timescale, set by the (3- 
decay half-lives. As a consequence the magic neutron num- 
bers (here N = 82) are reached at different astrophysical 
conditions and hence at different proton numbers, which is 
reflected in t he sh ifts of the abundance peaks (from Borzov 



and Goriely, 200C) 
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neutron-rich nuclei. Consequently the abun dance peaks 
are shifted to smaller A- values (see figure |2(| Borzov and 
Cor cry, [2000|). Si milar studies have been presented by 



Kratz et al (|998) 



matter by charge-current reactions, setting the proton- 
to-neutron ratio n/p or equivalently the Y e value of the 
r-process matter. A s shown by Fuller and Qian one ha s 
the simple relation (Qian, 1997; |Qian and Fuller, 1995) 



2. The possible role of neutrinos in the r-process 



Among the various possible astrophysical sites for the 



r-proc ess , the neutrino-d r iven w ind model (Witti et 



1994a, |; Wooslcy et al. , 1994 ) has attracted most at 



tention in recent years. Here it is assumed that the r- 
process occurs in the layers heated by neutrino emission 
and evaporating from the hot protoncutron star after 



core collapse in a type II supernova (Thompson et 
|200l| ). Adoptin g the parame ters of a supernova Simula 



tion by Wilson (Wilson, 1985) Woosley and collaborators 
obtained quite satisfying agreement between an r-process 
simulation and observation (Woosley et al, 1994). In 



the classical picture the r-process nuclides are made by 
successive capture of neutrons, starting from a seed nu- 
cleus with mass number A see d- Thus, to make the third 
r-process peak around A ~ 200 requires a large neutron- 
to-seed ratio of n/s ~ 200 — A scc d. In the Wilson super- 
nova models ( Wilsonj , 1985, 2001) this is achieved due to 
a high entropy found in the neutrino wind at late times 
(a few seconds after the bounce). However, other models 
with a different equation of state and treatment of dif- 
fusion do not obtain such high entropies; in these mod- 
els the r-process fails to make the A = 200 peak. To 
explain the strong sensitivity of the r-process nucleosyn- 
thesis on the entropy of the environment Qian (1997) 
noted that the slowest reaction in the nuclear network, 
which transforms protons, neutrons and a particles into 
r-process seed nuclei, is the 3-body a + a + n 
reaction. Due to its low binding energy (Eb 
MeV), 9 Be can be easily destroyed in a hot thermal en- 
vironment and thus the matter flow to nuclei heavier 
than 9 Be depends strongly on the entropy of the sur- 
rounding. The larger the entropy, the smaller the abun- 
dance of surviving 9 Be nuclei, which are then transformed 
into seed nuclei, and the larger the neutron-to-seed ra- 
tio. IVleyer (1995) pointed out that in a very strong neu- 
trino flux the slow 3-body a + a + n — > 9 Be reaction can 
be potentially bypassed by a sequence of two-body re- 
actions started by the neutrino-induced spallation of an 
a particle; e.g. 4 He(^, z/p) 3 H( 4 He,7) 7 Li( 4 He,7) n B an d 
4 He(j/,^n) 3 He( 4 He,7) 7 Be( 4 He,7) n C. This would speed 
up the mass flow to the seed nuclei and thus reduce the 
neutron-to-seed ratio. 



9 Be 
1.57 



Systematic stud i es by Hoffman et al. ( 1997 ) and 



Freiburghaus et al. (1999a) have shown that a successful 
r-process requires either large entropies at the Y e values 
currently obtained in supernova models, or smaller values 
for Y e . 

In the neutrino-driven wind model the extreme flux 
of v e and v e neutrinos from the protoneutron star inter- 
acts with the free protons and neutrons in the shocked 



p ~ L Ve (E Ve ) 



(9) 



As the neutrino energy luminosities are about equal 
for all species (L v ~ 10 52 erg s _1 ) the n/p-ratio is set 
by the ratio of average energies for the antineutrino and 
neutrino. As discussed above, their different opacities 
in the protoneutron star ensure that (Ep e ) > {E Vc ) and 
the matter is neutron-rich, as is required for a success- 
ful r-process. When the matter reaches cooler temper- 
atures, nucleosynthesis starts and the free protons are, 
in the first step, assemblied into a-particles, with some 
extra neutrons remaining. If these neutrons are still ex- 
posed to a large neutrino flux, it will change some of 
the neutrons into protons, which will then, together with 
additional neutrons, also be bound v ery quickly into a - 
particles. Thus, this so-called a effect ( Meyer et aL , 199S ) 
would severely reduce the final neutron-to-seed ratio and 
is therefore very counter-productive to a successful r- 
process. As mentioned above, the neutrino-nucleon cross 
sections are only considered to lowest order in super- 
nova simulations. The correction, introduced by the weak 
magnetism, acts to reduce the neutron-to-p r oton ratio in 
the neutrino-driven wind (Horowitz and Li, 1999|) . 

There are possible ways out of this dilemma: One so- 
lution is to remove the matter very quickly from the neu- 
tron star in order to reduce the neutrino fluxes for the 
a effect. Such a short dynamical timescale of the mate- 
rial in the wind is found in neutrino-driven wind models 
studied by Kajino and collaborators. These authors also 
observe that relativistic effects as well as nuclear reaction 
paths through neutron-rich lig ht elements m i ght b e help- 



ful for a successfu l r- process (Otsuki et 



sawa et 



triguing cure is discussed by 



2000 



Tera- 



2001a| Jb] IWanaio et al{ |200l|). A nother in 



McLaughlin et al. (1999) and 



Caldwell et al. ( |2000[ ) invoking matter-enhanced active- 
sterile neutrino oscillations to remove the v e from the 
r-process site. 

A simple estimate for the duration of the r-process can 
be had by adding up the half-lives of the waiting point 
nuclei, which results in about 1-2 seconds. However, 
in the neutrino-driven wind model it appears that the 
ejected matter passes through the region with the condi- 
tions suited for an r-process in shorter times (~ 0.5 s), 
implying that there might not be enough time for suffi- 
cient matter-flow from the seed to nuclides in the A ~ 200 
mass region. Such a 'time problem' is avoided if, as indi- 
cated above, the half-lives of the waiting point nuclei are 
shorter than conventionally assumed, or if, in a dynam- 
ically changing environment, the matter, which freezes 
out making the A ~ 200 r-process peak, breaks through 
the N = 50 and 82 waiting points closer to the neutron 
dripline, i.e. through nuclei with shorter half-lives, than 
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the matter which freezes out at these lower magic neutron 
numbers. 
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FIG. 27 Half-Lives of r-process waiting point nuclei with 
N = 50 (top panel), 82 (middle panel) and 126 (bottom panel) 
against charged-current (v e ,e~) reactions. For the neutrinos 
a Fermi-Dirac distribution with T = 4 MeV and zero chemi- 
cal potential (circles) and a luminosity of (L„ ~ 10 52 erg s _1 ) 
has been adopted. It is assumed that the reactions occur at 
a radius of 100 km, measured from the center of the neutron 
star. The half-lives can be significantly shorter if v e i/ MlT 
oscillations occur. The squares show the half-lives for neutri- 
nos with a Fermi-Dirac distribution with T = 8 MeV and zero 
chemical potential, which corresponds to complete v e f± f MlT 
oscillations. 

In an environment with large neutrino fluxes the 
matter-flow to heavier nuclei can also be sped up 
by cha rged- c urrent 



reactions (Nadyozhin and 
Panpv|, |1993 ; Qian et ai, 1997) which can compete with 
/3-decays. This is particularly important at the waiting 
point nuclei associated with N — 5 0, 82, and 126 . Fig- 



ure [27| shows the {ve,e ) half- lives ( Hektor et ai , 2000 



Langanke and Kolbe, 2001) for these waiting point nuclei, 
considering reasonable supernova neutrino parametriza- 
tions and assuming that the ejected matter has reached 
a radius of 100 km. Due to the dependence on L v the 
(v e ,e~) half- lives scale with r 2 . Indeed, a comparison 
with the (3 half- lives (figure ^) shows that {v ei e~) re- 
actions can be faster than the longest /3-decays of the 



N = 50, 82, and 126 waiting point nuclei, thus speed- 
ing up the breakthrough of the matter at the waiting 
points, if the r-process occurs rather close to the sur- 
face of the neutron star. This is even further enforced 
if v e u^ T oscillations occur due to the higher en- 
ergy spectrum of the supernova f MjT neutrinos. How- 
ever, the presence of charged-current reactions on nuclei 
in the neutrino-driven wind model implies also neutrino 
reactions on nucleons strengthening the a effect (Meyer 
" 1998| ). 




(Z+1.A) 

FIG. 28 Schematic view of the (v e ,e~) reaction on r-process 
nuclei. Due to the high neutrino energies the cross sections 
are dominated by transitions to the Fermi (IAS) and Gamow- 
Teller resonances. 

Under a strong neutrino flux the weak flow is deter- 
mined by an effective weak rate given by the addition of 
the charged-current (y e , e~) and the nuclear beta-decay 
rates ( [McLaughlin and Fuller . 1997] ). It has been ar- 
gued that the solar system r-process abundances pro- 
vide evidence for the weak steady-flow approximation, 
which implies that the observed abundances should be 
proportional to the ha lf-lives of th e ir progenit or nuclei 



on the r-process path (Kratz et al 



1993 



1988). If this 



is so, the r-process freeze-out must occur at conditions 
(i.e. radii) at which /3-decay dominates over (v ei e~) re- 
actions, at least for late times. The reason is that neu- 
trino capture on magic nuclei with N = 50, 82, 126 is not 
reduced if compared to the neighboring nuclides, as the 
capture occurs from a reservoir of neutrinos with suffi- 
ciently high energies to allow for transitions to the IAS 
and GT resonant states (see figure ^8|). Consequently 
(i/ e ,e~) cross sections scale approximately like the neu- 
tron excess (N — Z) , reflecting the Fermi and Ikeda sum- 
rules (see figure |27| ). However, the observed solar abun- 
dances around the A = 130 (N = 82) and 195 (N = 126) 
peaks do not show such a smooth dependence with A 
as will be the case if the effective weak rate is domi- 
nated by neutrino reactions. If we require that the de- 
cay half- lives dominate over the (y e , e~) reactions we can 
put constraints on the neutrino fluence in the neutrino- 
driven wind scenario, which is particularly strict if neu- 
trino oscillations occur. If, as an illustrative example, we 



31 



apply the constraint to the heavy waiting point nuclei 
with N — 126 (e.g. the nuclei with A ~ 199) and adopt 
the neutrino and beta half-lives from figure ^7|, /3-decay 
is only faster if the neutrino-nucleus reactions occur at 
distances larger than ~ 500 km. 



(figure |29|) . This shift depends on the neutrino exposures 
and allows constraints to be put on the total neutrino flu- 



ence i n the neutr i no-dr iven wind model (Haxton et 



1997 




Qian et al. 


, 1997) 
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The limits obtained this way 
are compatible with the values predicted by supernova 
models. Whether /3-delayed neutron emission, which has 
been neglected in (Haxton et al., 1997 ) mi g ht affect th e 
post-processing is an open question ( [Kratz , 2001 , 2002 ) . 

Attempts to include neutrino-induced reactions in the 
r-process network wit hin the neutrino-driven wind model 
have been rep orted inj Fuller and Mcyci L 1995| ; McLaugb 
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FIG. 29 Effect of post-processing by neutrino-induced reac- 
tions on the r-process abundance. The unprocessed distribu- 
tion (solid line) is compared with the distribution after post- 
processing (dashed line). A constant fiuence of T = 0.015 
has been assumed, which provides a best fit to the observed 
abundances for A = 183 — 87 (see inset). The observed abun- 
dance s are plotted as filled circles with error bars (from CJian 
et 4J |1997| )). 



As the Qp values in the very neutron-rich r-process 
nuclei are large, the IAS and GT resonant states reside 
at rather high excitation energies (~ 20-30 MeV) , in the 
daughter nuclei for (v e ,e~) reactions. This fact, com- 
bined with the small neutron separation energies in these 
nuclei, ensures that (v e ,e~) reactions, and also neutral- 
current (is, v') reactions, spa ll neutrons o u t of t he target 
nuclei ( [Haxton et al. , 1997| ; |Qian et~al, , 1997 ) (~ 5-7 
neutro ns fo r nuclei in the A = 1 95 mass region (Haxton 
et af], |1997 ; Hcktor et al, 2000| )). During the r-process, 
i.e. as long as the neutron source is strong enough to es- 
tablish (n, 7) (7, n) equilibrium, the neutrons will im- 
mediately be recaptured leading to no effect on the abun- 
dance distribution. However, once the neutron source has 
ceased, e.g. after freeze-out, and if the r-process matter in 
the neutrino-driven wind model is still subject to strong 
neutrino fluxes, neutrons liberated during this phase by 
neutrino-induced reac tions, will not be re c aptured and 



this post-processing (Haxton et al., 1997; Qian et al. 



1997) leads to changes in the r-proc e ss ab undance dis- 



tribution. It is argued ( Haxton et al , 1997 ) that due to 
the smooth dependence of the neutrino cross sections on 
the mass number, the post-processing, in general, shifts 
abundances from the peaks to the wings at lower A- values 



lin and Fuller, 1995; Tcrasawa et 



Mcycr|, |l9bta| ; [Mcbaugn- 
q/J,' |2001a| JbP- In partic- 



ular Meyer et al. (1998) studied the competition of the 
a-effect with the possible speed-up of the matter-flow by 
charged-current reactions on nuclei. These authors esti- 
mated the respective charged-current cross sections for 
allowed transitions on the basis of the independent par- 
ticle model. Improving on this treatment, Borzov and 
Goriely calculated (v e ,e~) cross sections for supernova 
neutrinos (Fermi-Dirac distribution with T — 4 MeV) 
within the ETFSI method, consiste ntly with their most 
recen t estimates for the (3 half-lives (Borzov and Goriely, 
200C ) . RPA-based neutrino- nucleus cross sections for se- 



lected nuclei have been re ported in (Hektor et al., 2000 
Burman and Engel , 1998 ). Very recently a tabulation 
with charged- and neutral-current total and partial neu- 
tron spallation cross sections hav e become available for 
the neutron -rich r-process nuclei (Langanke and Kolbc, 
200 1| , [2002D . This tabulation is based on the RPA and 
considers allowed and forbidden transitions. Further- 
more, the cross sections are tabulated for various super- 
nova neutrino distributions, thus also allowing study of 
the influence of complete neutrino oscillations on the r- 
process. 



Terasawa et al. ( 2001a ^) have performed studies simi- 
lar to the pioneering work of Meyer et al. (1998), however, 



using their neutrino-driven wind model and th e complete 



and Kolbe 




2001 


, |2002) 



A typical result is shown in 
figure |30|, where, however, the simplifying, but reason- 
able assumption of an exponential time-dependence of 
the matter- flow, governed by the parameter Tdyn, away 
from the neutron star has been assumed. The quan- 
tity (A) defines the average mass of the heavy seed nu- 
clei present at the beginning of the r-process, defined at 
T = 2.5 x 10 9 K. The neutron-to-seed ratio n/s is very 
sensitive to the dynamical evolution time. This comes 
about as the shorter Td yn , the less time is available to 
assemble the seed nuclei from a-particles and neutrons. 
Consequently the abundance of seed nuclei decreases for 
shorter Td yn , increasing the n/s ratio. If the neutrino 
flux is artificially switched off, matter-flow to the 3rd r- 
process peak at A ~ 200 (2nd r-process peak at A ~ 130) 
is achieved if Td yn < 0.01 s (rd yn < 1 s). The consistent 
inclusion of neutrino reactions is counter-productive to a 
successful r-process. This effect becomes more dramatic 
if the matter-flow is slow as then the a-effect strongly 
suppresses the availability of free neutrons at the begin- 
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FIG. 30 The lower panel shows the average mass number of 
heavy seed nuclei (A) and the neutron-to-seed ratio (n/Seed) 
in a neutrino-driven wind simulation with an exponential time 
dependence of the matter flow, determined by the parameter 
Tdyn- The open squares show the results for a simulation 
including neutrino-nucleus reactions, while the full triangles 
refer to a study in which the neutrino-reactions on nuclei have 
been switched off. The upper panel shows the sum of the mass 
number of the seed nucleus plus the neutron-to-seed ratio. 
This quantity shows up to which mass number the r-process 
can produce nuclides. The horizontal lines indicate the second 
(A ~ 130) and third (^4 ~ 200) r-process peaks as well as the 
position of the smaller r-process peak related to the deformed 
nuclei in the rare-earth region (REE peak), (courtesy of M. 
Terasawa) . 



ning of the r-process (see also Meyer et al. , 1998| ). No 
r-process, i.e. no production of nuclides in the second r- 
process peak at A ~ 130, is observed if Tdy n >0.05 s. 



Recently Qian ( 2002 ) has suggested that, within the 
neutrino-driven wind r-process scenario, charged-current 
neutrino reactions can induce fission reactions on r- 
process progenitor nuclei heavier than lead and that the 
fission products account for the observed r-process abun- 
dance in metal-poor stars (Bncdcn et al, 2000). These 



observed abundances show a peak around mass number 
A ~ 195, which follows the solar r-process distribution, 
and enhanced structures at aro und A ~ 90 and ~ 132 
which, due to the suggestion of Qian ( [2002 ), are fission 
products, which do not follow a solar r-process pattern. 
First calculations of neutrino-i nduced fission cro ss sec- 
tions have been performed by (Kolbe et al., 2002| ) using 
a combination of RPA model, to calculate the (i/ e ,e _ ) 
excitation function, and statistical model to determine 
the final branching probabilities. The neutrino-induced 
fission cross sections (see figure [u]) are quite large as 
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FIG. 31 Total (v e ,e~) (circles) and partial (v e ,e~n) (trian- 
gles) and neutrino-induced fission cross sections (squares) for 
selected thorium and uranium isotopes. The calculations have 
been performed for Fermi-Dirac v e distributions with tem- 
perature T = 4 MeV and 8 MeV. The first reflects a typical 
supernova v e spectrum, while the latter assumes complete 
Va t — » neutrino oscillations. 



the progenitor nuclei are neutronrich, which increases 
the Fermi and Gamow- Teller contributions to the total 
cross sections and places their strengths at energies above 
the fission barrier in the daughter nucleus. The calcula- 
tions shown in figure pl|use the fission barriers derived by 
980|). Modern evaluations predict 
200 ID which 



Howard and Moller 



larger fission barriers (e.g. Mamdouh et 



would reduce the fission cross section. The cross sections 
can be significantly enlarged if z^. T — > v e neutrino oscil- 
lations occur during the neutrino-driven wind r-process 
scenario. 



VII. THE NEUTRINO PROCESS 

When the flux of neutrinos generated by the cooling 
of the neutron star passes through the overlying shells of 
heavy elements interesting nuclear transmutations are in- 
duced, despite the small neutrino-nucleus cross sections. 
Of particular interest here are neutral-current reactions 
as they can be induced by v^, v T neutrinos and their an- 
tiparticles with the higher energy spectra {{E„) ~ 25 
MeV) . These neutrinos are energetic enough to excite the 
GT resonant state, and more importantly, also the giant 
dipole resonant states. These states usually reside above 
particle thresholds and hence decay mainly by proton or 
neutron emission, generating new nuclides. The neutrino 
reaction rates are too small to affect the abundances of 
the parent nuclei, but they can noticeably contribute to 
the production of the (sometimes much less abundant) 
daughter nuclei. As a rule-of-thumb, the neutrino pro- 
cess, i.e. the synthesis of nuclides by neutrino-induced 
spallation in a supernova, can become a significant pro- 
duction process for the daughter nuclide if one wants to 
explain abundan ce ratios of pare nt-to-daughter which ex- 
ceed about 10 3 dWoosleyj , |200l|) . 

The most interesting neutrino nucleosynthesis occurs 
in the outer burning shells of the massive star which have 
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not been affected by the fatal core collapse in the cen- 
ter when the neutrinos pass through. However, a little 
later the shock reaches these shells and the matter will be 
subjected to rather high temperatures initiating fast nu- 
clear reactions involving also the nuclides just produced 
by the neutrino-induced reactions. Hence studies of the 
neutrino process depend on various neutrino-nucleus re- 
action rates and on the neutrino spectra and fluxes, es- 
pecially of the v^, v T neutrinos, and they require a rather 
moderate nuclear network to simulate the effects of the 
re-processing by the shock. The first inv estigation of the 
neutrino process has been reported in ( Woosley et al , 
1990| ). A more recent study ([Woosley and Weaver|, |1995|) 



confirmed the main result that a few specific nuclei ('Li, 
n B, 19 F) are being made by the neutrino process in sig- 
nificant fractions. For example, n B and 19 F are being 
made by [y, v'p) and (y, v'n) (followed by a /3-decay) re- 
actions o n the abundant 1 2 C an d 20 Ne, respectively. As 
noted by Woosley et al. ( |L990 ) neutrino nucleosynthe- 
sis can also contribute to the production of 180 Ta (see 
above) and 138 La. First calculations of the relevant to- 
tal and parti al neutrino- in d uced cross sectio ns have been 
reported in (Heger et al, 2002), (see also Belie et al, 
2002 ) . The nucleus 138 La is special as it is mainly made 
by the charged-current reaction 138 Ba(^ e , e~) 138 La, while 
the lighter nuclei ( 7 Li, 11 B, 19 F) and dominantly also 
180 Ta, are being produced by neutral-current reactions 
induced by ^,^ r neutrinos. The neutrino process is 
therefore sensitive to v e and (y^, v T ) neutrinos, which are 
expected to have different spectra in Type II supernovae 
(see section V.C), and hence it can test this prediction. 
Moreover, the 138 La nucleosynthesis is sensitive to neu- 
trino oscillations, as this nuclide would be significantly 
overproduced if supernova , v T neutrinos have a no- 
ticeably larger average energy and if they oscillate into 
v e neutrinos before reaching the 138 La production site 
(helium shell) in massive stars. 

Neutrino-induced nucleon spallation on 12 C can also 
knock-out a deuteron or a proton-neutron pair, in this 
way producing 10 B. The expected 10 B/ 11 B abundance 
ratio in neutrino-nucleosynthesis is ~ 0.05, which is sig- 
nificantly smaller than the observed abundance ratio, 
0.25 ( Haxton| , 2000). Thus there must be a second pro- 
cess which contributes to the production of 10 ' n B. These 
are reactions of energetic protons on 12 C in cosmic rays, 
which yield a ratio of 10 B/ 11 B of about 0.5, larger than 
the observed value. A solution might be that the two nu- 
clides are being produced by both mechanisms, neutrino- 
nucleosynthesis and cosmic ray spallation. It is interest- 
ing to note that the first process, being associated with 
supernovae, is a primary process, while the latter is a 
secondary process, as it requires the existence of protons 
and 12 C in the interstellar medium. As a consequence the 
10 B/ 11 B abundance ratio should have changed during the 
history of the galaxy. This can be tested once observers 
are able to distinguish between the abundanc es of the 
two different boron nuclides in stellar spectra (Haxtor, 
2000| ). 



VIII. BINARY SYSTEMS 

Weak processes can also play interesting roles in the 
evolution and nucleosynthesis processes in close binary 
systems where one component is a compact object (white 
dwarf or neutron star) and the other a massive star. If 
the latter expands during helium core burning, mass flow 
from the hydrogen envelopment of the star onto the sur- 
face of the compact object sets in. If the respective mass 
accretion rate is rather low (10 _8 -10 -10 M Q y^ 1 ), the 
hydrogen, accreted on the surface of the compact object, 
burns explosively under degenerate conditions leading to 
a nova (if the compact object is a white dwarf) or an 
x-ray burst (neutron star). This means that the energy 
released by the nuclear reactions is used to heat the mat- 
ter rather than for expansion. The rise in temperature 
increases the nuclear reaction rates, triggering a ther- 
monuclear runaway until the degeneracy is finally lifted 
and an outer layer of matter is ejected. In a type la 
supernova the faster mass accretion rate on to the sur- 
face of a white dwarf (likely a carbon-oxygen white dwarf 
with sub-Chandrasekhar mass ^0.7 M©) leads to steady 
hydrogen and subsequently helium burning. If the grow- 
ing mass of the white dwarf exceeds the Chandrasekhar 
mass, contraction sets in and the carbon in the center 
ignites by fusion reactions with screening enhancement. 
As the environment is highly degenerate, a thermonu- 
clear runaway is triggered which eventually will explode 
the entire star. 



A. Novae 

The main energy source of a nova is the CNO cycle, 
with additional burning from proton reactions on nuclei 
between neon and sulfur, if the white dwarf also con- 
tained some 20 Ne (Truran, 1983). A nova expels matter 
which is enriched in (3- unstable i4 ' 15 for carbon-oxygen 
white dwarfs (leading to the production of stable 14,15 N 
nuclides) or can contain nuclides up to the sulfur mass 



region for neon - oxygen nov ae (lose and Hernanz, 1998 



(Starrfield et al , 2000 , 1998 ). The most important role of 



weak-interaction processes in novae is their limitation of 
the energy generation during the thermonuclear runaway. 
An interesting branching occurs at T ~ 8 x 10 7 K. Fol- 
lower temperatures the /3-decay of 13 N with a half-life of 
~ 10 m dominates over the 13 N(p, 7) reaction and sets the 
timescale for the nuclear burning. As charged-particle 
fusion reactions are sensitively dependent on tempera- 
ture, their reaction rates strongly increase with rising 
temperature and for T > 8 x 10 7 K and densities of or- 
der 10 4 g cm~ 3 the proton capture on 13 N is faster than 
the /3-decay. The CNO cycle turns into the hot CNO 
cycle and now the positron decay of 15 O with a half- 
life of 122 s is the slowest reaction occuring in the CNO 
nova network. It turns out that, once the degeneracy is 
lifted, the dynamical expansion timescale is faster than 
the one for nuclear energy generation, set by the 15 O half- 
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life. As a consequence the runaway is quenched ( Truran 
1982 ) . We mention that the determination of the domi- 



nant resonant contribution to the 13 N(p, 7) 14 rate has 
been the first successful application of the Coulomb dis- 
sociat ion te chnique in nuclear astrophysics ( Motobayashi 
et ofl, |1991| ). 

Weak-interaction rates relevant for nova studies can 
be derived from either experimental data or shell-model 
calculations. 



B. X-ray Bursts 

An x-ray burst is explained as a thermonuclear run- 
away in th e hydrogen-r i ch envelope of an a c cretin g neu- 



tron star (Lewin et ai, 1993; Faam et ai, 1993). Due 



to the higher gravitational potential of a neutron star, 
the accreted matter on the surface reaches larger densi- 
ties than in a nova (up to a few 10 6 g cm~ 3 ) and the 
t emperature during the runaway can rise up to 2 x 10 9 K 



( [Schatz et ai , 1998] ). Under these conditions hydrogen 



burning is explosively fast. The trigger reactions of the 
runaway are the triple-alpha reactions and the break- 
out from the hot CNO cycle, mainly by a-capture on 
15 and 18 Ne. We note that these two reaction rates 
are yet insufficiently known due to uncerta i n res o nant 
contri bution at low energies flGorrcs et~al\ , [1995| ; [Mao] 
et af], [1996 ). The thermonuclear runaway is driven first 
by the ap process, a sequence of (a,p) and (p, 7) reac- 
tions which shifts the ashes of the hot CNO cycle to the 
Ar and Ca mass region, and then by the rp-process (short 
for rapid proton capture process). The rp-process repre- 
sents proton-capture reactions along the proton dripline 
and subsequent /?-decays of dripline nuclei processing the 
mate rial from the Ar-Ca region up to 56 Ni flSchatz et al 



1998). The (i half-lives on the rp-process path up to ^Ni 



are known sufficiently well. 

The runaway freezes out in thermal equilibrium at 
peak temperatures of (2-3) x 10 9 K with an abundance 
distribution rich in 56 Ni, forming Ni oceans on the sur- 
face of the neutron star. Further matter flow at these 
temperatures is suppressed due to the low proton sep- 
aration energy of 57 Cu and the long half-life of 56 Ni. 
Re-ignition of the rp-process then takes place during the 
cooling phase, starting with proton capture on 56 Ni and 
potentially shifting matter up to the 100 Sn region where 
the rp -proc ess ends in a cycle in the Sn-Te-I range ( Bchatz 
et a[ , 2001 ) . A matter flow to even heavier nuclei is possi- 
ble, if the rp-process operates in a repetitive mannor; i.e. 
a new rp-process is ignited after the ashes of the previous 
process have decayed to stability and before these nuclei 
have sunk too deep into the crust of the neutron star (see 
below) . Such repetitive rp-process m odels, shortly calle d 
(rp) 2 -process, have been studied by ( Boyd et al , [2002 ). 

The reaction path beyond 56 Ni runs through the even- 
even N = Z nuclei which due to their known long half- 
lives ( 64 Ge has a half-life of 63.7(25) s) represent a strong 
impedance to the matter-flow. This can also not be over- 



come always by proton captures as for some of the a- 
like nuclei the resulting odd-yl nucleus is proton-unbound 
and exists only as a resonance. Such situations occur, 
for example, for the 68 Se(p, 7) 69 Br and 72 Kr(p, 7 ) 73 Rb 
reactions. It has been suggested (Gorres et al, 199?;) 



that the gap in the reaction path can be bridged by two- 
proton captures, with the resonance serving as interme- 
diate state (like in the triple-a reaction). The reaction 
rate for such a two-step process depends crucially on the 
resonance energy, with some appropriate screening cor- 
rections. 

Most half-lives along the rp-process path up to the 
80 Zr region are known experimentally. This region of the 
nuclear mass chart is known for strong ground state de- 
formations, caused by coupling of the p/-shell orbitals to 
the gg/2 and d§/2 levels. The strong deformation makes 
theoretical half-life predictions quite inaccurate, mainly 
due to uncertainties in the Q-valucs stemming from in- 
sufficiently wellknown mass differences. We note that 
the effective half-life of a nucleus along the rp-process 
path could be affected by the feeding of isomeric states 
in the proton captures or by thermal population of ex- 
cited states in general. Again, deformation plays a ma- 
jor role as then even in a-like nuclei excited states are 
at rather low energies (e.g. the first 2 + state in 80 Zr is 
at 290 keV). Apparently a measurement of the half- lives 
is indispensable for rp-process studies beyond A = 80. 
An important step has recently been taken by measuring 
t he half-life of 80 Zr at the Holifield Facility in Oak Ridge 
(Ressler et al , 2000] ). The experimental value of 4.lto'ej s 
reduces the previous (theoretical) uncertainty consider- 
ably and, in fact, it is shorter than the value adopted pre- 
viously in x-ray bursts simulations. Fast proton captures 
on the daughter products 80 Y and 80 Sr allow matter- flow 
to heavier nuclei, with the a-nucleus 84 Mo (N = Z = 42) 
being the next bottleneck. Experiments to measure this 
important half-life are in progress. 

The nucleosynthesis during the cooling phase in an 
x-ray burst alters considerably the abundance distribu- 
tion in the atmosphere, ocean and crust of the neutron 
star. For example, the rp-process may be a possible con- 
tributor to the presently unexplained relatively high ob- 
served abundance o f light p- nuclei like 92 Mo and 96 Ru 
flBchatz et ai , 1998 ). This assumes that the matter pro- 



duced in the x-ray burst gets expelled out of the large 
gravitational potential of the neutron star, which is still 
questionable. Due to continuing accretion the rp-process 
ashes are pressed into the ocean and crust of the neu- 
tron star, replacing there the neutron star's original ma- 
terial. When the ashes sink into the crust, they reach 
regions of higher densities, and relatedly, larger electron 
chemical potentials. Thus, consecutive electron captures 
will become energetically favorable and make the ashes 
more neutron-rich. At densities beyond neutron drip 
(p ~ (4-6) x 10 11 g cm~ 3 ) neutron emissions become 
possible and at eve n higher densities pycnonuclear re- 
actions can set in ( Bisnovatyi-Kogan and Chechetkm, 



1979] ; |Haensel and Zdunik| , |1990| ; |Satc| , [19791 ). Impor- 
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tantly, these processes (electron capture, pycnonuclear 
reactions) generate energies which can be locally stored 
in the neutron star' s ocean and crust and will affect their 
thermal properties ( Brown and Bildstcij , (19981) . .Previous 
studies of these processes in accreting neutron stars have 
assumed that iron is the endproduct of nuclear burning 
and the s ole nucleus reaching t he cr ust of the neutron 
star (e.g. Haensel and Zdunik , 1990 ). But clearly the 
r p-process produces a wide mixture of heavy elements 



( [Schatz et al. , 1999), where the abundance distribution 



depends on the accretion rate. 

The ashes consist mainly of even-even N = Z nuclei, 
for which electron capture at neutron star conditions oc- 
cur always in steps of two. At first, the capture on the 
even-even nucleus sets in once sufficiently high-energy 
electrons are available to effectively overcome the Qec~ 
value to the odd-odd daughter nucleus. Due to nuclear 
pairing, which favors even-even nuclei, the QEC-vahie for 
the produced daughter nucleus is noticeably lower so that 
electron capture on the daughter readily follows at the 
same conditions. The energy gain of the double-electron 
capture is of order the difference of the two QEC-v&hies; 
this gain is split between the emitted neutrino and a lo- 
cal heating. Considering a blob of accreted matter ini- 
tially consisting solely of 56 Ni and assuming temperature 
T = 0, the evolution of this blob was followed on the neu- 
t ron star surface unti l neut ron-drip densities and beyond 
( Haensel and Zdunik , 199C). The rp-process simulations, 
however, indicate a finite temperature of the ashes of a 
few 10 s K, allowing electron capture already from the 
high-energy tail of the electron distribution and signifi- 
cantly reducing the required densities. 



C. Type la Supernovae 

Type la supernovae at high redshifts serve currently 
as the standard candles for the largest distances in the 
universe. Importantly recent surveys of such distant su- 
pernovae provide evidence for an acceleratin g expansion 
of the universe over the last several 10 9 years (Pcrlmuttcr 
et a{], |1999| ; |Riess et 



1998) 



Type la supernovae have been identified as thermonu- 
clear explosions of accreting white dwarfs with high ac- 
cretion rates in a close binary system. While the general 
explosion mechanism is probably understood, several is- 
sues are still open like the masses of the stars in the bi- 
nary or the carbon/oxygen ratio and distribution in the 
white dwarf. The probably most important problem yet 
is the modelling of the matter transport during the explo- 
sion and the velocity of the burning fr ont, both requir 
ing multidimen s ional simulations (e.g 



Niei freyer| , [2000| ; |Reinecke et al\ , |1999| ; |Woosley| , |200l|) 
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FIG. 32 F e profile as a function of rad ial mass for th e stan- 
dard type la supernova model WS15 ( (Sfomoto et al. , 1984) 
using the FFN and the shell-model weak-interaction rates 
(LMP), (courtesy of F. Brachwitz). 



tor of two compared with the relative solar abundances. 
This requirement puts stringent constraints on models, 
in particular on the central den sity of the progenito r 
white dwarf and the flame speed (Iwamoto et al, 1999). 
When the flame travels outwards, it heats the matter 
to temperatures of a few 10 9 K and brings its compo- 
sition close to nuclear statistical equilibrium (NSE). As 
the original matter ( 16 0, 12 C) had an electron-to-baryon 
ratio of Y e = 0.5, the NSE composition is dominated by 
56 Ni, which after being expelled decays to 56 Fe. However, 
behind the flame front, which travels with a few per- 
cent of the local sound speed ( Niemeyer and Hillebrandt , 



1995), electron captures occur, which lower Y e and drive 
the matter composition more neutron-rich. This effect 
is larger the greater the central density of the white 
dwarf (which increases the electron capture rates) and 
the slower the flame speed (which allows more time for 
electron captures) . Figure shows the Y e profile ob- 



taine d in a standard type la model WS15 (Nomoto et 



1984), with slow deflagration flame speed (1.5% of sound 



velocity), central ignition density p = 2.1 x 10 9 g cm"' 
and a transition from deflagration to detonation at den- 
sity p = 2. 1x 10 7 g cm' 3 . Th e calculations h ave been per- 



formed by Brachwitz (2001) with the FFN (Fuller et 



1982b) and shell model (Langanke and Martmez-Pincdo 



Type la supernovae contribute about half to the abun- 
dance of Fe-group nuclides in galactic evolution. Thus 
one can expect that type la supernovae should not over- 
produce abundances of nuclides in the iron group, such model, it increases the central density by about a fac- 
as 54 Cr or 50 Ti, relative to 56 Fe by more than a fac- tor of 1.3 (iBrachwitz et all 120001; IWoosleyj, 120011). This 



20011) weak interaction rate sets. The differences are 
quite significant, even if one considers that about 60% of 
the captures occur on free protons, which are unaffected 
by the differences in these rate sets. Under otherwise 
identical conditions the slower shell model rates yield a 
central Y e value of 0.45, which is about 0.01 larger than 
for the FFN rates. Consequently very neutron-rich nu- 
clei with Y e < 0.45 are significantly suppressed (see fig- 
ure^). In fact, no nuclide is significantly overproduced 
in this model compared to the solar abundance. The net 
effect of the new rates is that, for an otherwise unchanged 
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FIG. 33 Ratio of calculated to solar abundances predicted by the WS15 model ( Nomoto et al. , 1984 ) using the FFN and the 
shell-model rates (LMP). The ordinate is normalized to 56 Fe. Intermediate-mass elements exist but are underproduced by a 
factor of 2-3 for SNe la models in comparison to Fe group elements. Using the FFN rates the Fe group does not show a 
composition close to solar. Especially 50 Ti and 54 Cr are strongly overproduced by more than a factor of 3. The change from 
FFN rates (left) to LMP (right) reduces the overproduction over solar close to the acceptable limit of a factor 2 (courtesy of 
F. Brachwitz). 



can have quite interesting consequences if one wants to 
use the nucleosynthesis constraint to distinguish between 
two quite distinct type la models. On the basis of recent 
models it has been concluded that the majority of type 
la progenitors grow towards the Chandrasekh ar mass 
through steady hydrogen and helium burning (Hachisu 
et aQ [1999 ). Such systems would lead to rather low cen- 
tral densities p < 2 x 10 9 g cm -3 . In these models, only 
a small fraction of progenitors would deviate from steady 
hydrogen burning at the end of the accretion history ex- 
periencing weak hydrogen flashes; such cases correspond 
to the W-model discussed above yielding higher central 
densities. We stress, however, that changes in the nucle- 
osynthesis caused by differences in the central densities 
in the models can be counterbalanced by changes of the 
flame speed. 



IX. CONCLUSIONS AND FUTURE PERSPECTIVES 

It has long been recognized that nuclear weak- 
interaction processes play essential roles in many astro- 
physical scenarios. In a few cases these are specific re- 
actions, which are particularly important, and these re- 
actions have then been studied with increasingly refined 
models. Examples are the various weak-interaction pro- 
cesses in the solar hydrogen bu rning chains including the 
initi al p + p fusion reaction ( Kamionkowski and Bah 



calf! |1994| ; |Kong and Rayndaffiobll ; park et al\, |1998[ 
2001a|) or the very challenging 3 Hc + p reaction (Marcucci 



et qjj, |2000| , |200lt |Park et al\ , |2001b]) , which generates the 
highest-energy neutrinos in the sun. Another typical ex- 
ample is the solar electron capture rate on 7 Be, where 
the nuclear matrix element can be determined from the 
experimental lifetime of atomic 7 Be, while the proper de- 
scription of the solar plasma effects on the capture rate 



with the desired accuracy has been quite demanding. 

However, most astrophysical applications require the 
knowledge of weak-interaction rates for a huge body of 
nuclei. If, like for s-process nucleosynthesis, the nuclei in- 
volved are close to the valley of stability and hence quite 
long-lived, the needed rates (usually half-lives) have been 
determined experimentally in decade-long efforts. For 
the s-process the challenge now focuses on the branching 
point nuclei where the reaction flow branches into two 
(or more) paths and where, in some cases, the observed 
relative abundances of nuclides along the different paths 
depend on the stellar conditions (temperature, neutron 
density) and hence allow determination of these quanti- 
ties inside the star. Again, the approach is to measure the 
necessary data, e.g. half-lives of excited nuclear states. 

Other astrophysical scenarios involve nuclei far-off sta- 
bility, often under extreme conditions (high density, neu- 
tron flux, temperature). These astrophysical sites in- 
clude core-collapse (type II) and thermonuclear (type 
la) supernovae, r-process nucleosynthesis and explosive 
hydrogen burning; and the interest in all of these has 
been boosted recently by novel observations and data 
(supernova 1987A, high-redshift supernova survey, Hub- 
ble Space Telescope, . . .). A direct experimental deter- 
mination of the respective stellar weak-interaction rates 
has been possible in a few cases, like the half-lives of r- 
process and rp-process waiting point nuclei. However, in 
nearly all cases the weak-interaction processes had to be 
theoretically modeled so far - a very demanding job if one 
considers that often results for many hundreds of nuclei 
for a large range of stellar conditions are needed. These 
data were then derived globally based on parametrized 
nuclear structure arguments, as an appropriate treatment 
of the involved nuclear structure problem was prohibited 
by both the available computational capabilities and the 
lack of experimental guidance. Although evaluation of 
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rate sets for astrophysical purposes often appears to be 
a theoretical problem, the second point - experimental 
guidance - is crucial and often overlooked. It is Willy 
Fowler's strategy and legacy that nuclear models used to 
derive nuclear ingredients in astrophysical applications, 
should be consistent, but more importantly they should 
be accurate and, as a consequence, experimental data 
are to be used whenever available. Therefore the role 
of experimental data in nuclear astrophysics is twofold: 
If possible, they supply the needed information directly, 
or equally important, they constitute constraints and 
guidance for the nuclear models from which the needed 
information is obtained. Thus, the renaissance of nu- 
clear structure, which we have witnessed in recent years, 
has two consequences in nuclear astrophysics. The re- 
cent development of new facilities, techniques and de- 
vices brought a large flood of new experimental informa- 
tion, in particular for the proton- and neutron-rich nuclei 
away from stability. These data indicate that the nuclear 
structure models adopted to derive the global data and 



tance to the source and the measured intensity of the 
7-line allows determination of how much 44 Ti has been 
ejected into the interstellar medium by the supernova 
event. Furthermore, as 7-rays can, in contrast to optical 
wave lengths, escape from the galactic bulk, 7-ray obser- 
vation allows to detect histori cal supcrnovac whic h have 
not been observed op tically (Tyudin et al , 1998[ ), 



(also 

observed in x-rays by Aschcnbach| , 199S ). Such searches 



for historical supernovae and an improved determination 
of the supernova frequency in our galaxy will be one of 
the main missions of future 7-ray observatories in space, 
like INTEGRAL 3 . A longer-lived radioactive nuclide pro- 
duced in supernovae is 60 Fe. Investigations of rock sam- 
ples taken from the ocean floor, by precision accelerator 
mass spectroscopy, found a significant increase of 60 Fe 
abundance compared to other iron isotopes p ointing to a 



close -by supernova about 5 million years ago ( Knie et al 



1999) 



rateisets were usually too simple and improvements were 
warranted. The experimental renaissance went hand-m- 
hand with decisive progress in nuclear structure theory, 
made possible by the development of new models and 
better computer hardware and software. Due to both 
experimental and theoretical advances, it becomes now 
possible to calculate nuclear data sets for astrophysical 
applications on the basis of realistic models rather than 
on crude and often oversimplified parametrizations. This 



Weak-processes on nuclei and electrons are the means 
to observe sola r and supernova neutrinos (Balantckin and 



Haxton , |199S|). Solar neutrinos have rather low energies 



(E v < 14 MeV) and hence induce specific low-lying tran- 
sitions which are theoretically modelled best by shell- 
model calculations. A pplications have be e n perf ormed 
for example, for 37 CI ( Aufderheidc et al , 1994a|) 



40 



Ar 



flOrmand et al\ , |l995| ) and yl Ga ( |Haxton|, |1998|), the 



detec tor material in the Homestake (Cleveland et al 
ICARUS 4 , and GALLEX/SAGE/GNO de t ectors 



199S 



(lATxlurashitov et al\ 1999 ; [Altmann et al\ , |2000| ; [Ham 



tions 

The advances in nuclear structure modelling have also 
lead to progress in astrophysically important nuclear in- 
put other than weak-interaction processes. Typical ex- 
amples are the equation of state, derived on the ba- 
sis of the relativistic mean-field model guided by the 



revi fw presents a summary ot recent tneoretical caicma- pel et g | 1999 | ) ; reipectively. Supernova neutrinos have 



relativ istic Bruckner-Hartree-Fock theory ( Shen et al 
1998a| H), which serves as an alternative to the standard 



higher energies; in particular, the energies of fj,, r neutri- 
nos and antineutrinos are expected to be high enough to 
excite the giant dipole resonances in nuclei. Studies for 
detector materials like Na, Fe, Pb have been performed 
in hybrid approaches combining shell model calculations 
for allowed transitions with the RPA model for forbid- 
den transitions or treating the forbidden transitions on 



Lattimer-Swesty EOS (Lattimcr and Swesty, 1991), or 
the nuclear mass table and level density parametriza- 
tions determined within the framework of the Hartree- 



the basis of the Goldhabe r - Telle r model (Fuller et al 



1999; Kolbe and Langankc , 2001) 



Fock model with BCS pa iring (Demetriou and Goriefy 



2001 ; Goriely et al , 2001 ). Such mass tables, or equiva- 
lently neutron separation energies, play an essential role 
in r-process nucleosynthesis and are conventionally de- 
rived by parametrizations constrained to known masses. 

Other astrophysical areas, loosely or indirectly related 
to the topic of this review, have also benefitted from the 
progress in modelling nuclear weak-interaction processes. 
A field with rapidly growing importance is Gamma-ray 
Astronomy with beta-unstable nuclei. Due to 7-ray ob- 
servatories in space it has been possible in recent years 
to search the sky for sources of known 7-rays which can 
then be associated with recent nucleosynthesis activities. 
A highlight has been the observation of the 1.157 MeV 7- 
line, produced in the /3-deca y scheme of 44 Ti, in the Cas- 

Know- 



Despite the experimental and theoretical progress, lack 
of knowledge of relevant or accurate weak-interaction 
data still constitutes a major obstacle in the simulation 
of some astrophysical scenarios today. This refers mainly 
to type II supernovae and r-process nucleosynthesis. 

Core-collapse supernovae, as the breeding places of 
carbon and oxygen, and hence life, in the universe, at- 
tract currently significant attention and the quest to 
definitely identify the explosion mechanism is on the 
agenda of several international and interdisciplinary col- 
laborations. While most of the efforts concentrate on 
computational developments towards a multidimensional 
treatment of the hydrodynamics and the neutrino trans- 
port, some relevant and potentially important nuclear 



siopeia supernova remnant (lyudin et al. 
ing the date of the supernova, the 44 Ti half-life, the dis- 



199J) 



3 http: / /astro. estec.esa.nl/SA-general/Projects/Integral/integral. html 

4 http:/ /www. aquila.infn.it/icarus 



38 



problems remain. The improved description of weak- 
interaction rates in the iron mass range has lead to sig- 
nificant changes in the presupernova models. Evaluation 
of electron capture rates for heavier nuclei will be avail- 
able in the near future. First results imply that capture 
on heavy nuclei, usually ignored in collapse simulations, 
can compete with the capture on free protons. Whether 
the inclusion of this process then leads to changes in 
the collapse trajectory has yet to be seen. Further, fi- 
nite temperature neutrino-nucleus rates will also soon be 
available. This will then test whether inelastic neutrino 
scattering on nuclei, yet not modelled in the simulations, 
influences the collapse dynamics or supports the revival 
of the shock wave by preheating of the infalling matter 
after the bounce. The largest nuclear uncertainty in col- 
lapse simulations are likely associated with the descrip- 
tion of nuclear matter at high density, extreme isospin 
and finite temperature. In particular a reliable descrip- 
tion of the neutrino opacity in nuclear matter might very 
well be what is needed for successful explosion simula- 
tions. Nucleon-nucleon correlations have been identified 
to strongly influence the neutrino opacities and nuclear 
models like the RPA are quite useful to guide the way. 
Ultimately one would like to see the many-body Monte 
Carlo techniques, which have been so successfully ap- 
plied to few-body systems or the nuclear shell model, ex- 
tended to the nuclear matter prob lem. First steps along 
t hese lines have been rep o rted b y Schmidt and Fantoni 
( |f 999| ) and Fantoni et al\ ( |2001 ), who proposed a novel 
constrained-path Diffusion Monte Carlo model to study 
nuclear matter at temperature T — 0. The Shell Model 
Monte Carlo model, formulated in momentum space and 
naturally at finite temperature, constitutes an alterna- 
tive approach. First attempts in this direction have been 



taken by the Caltech group (Zheng, 1996) and by Rom ■ 
bou||| ( [1998 ). Whether the notorious sign-problem in the 



SMMC approach can also be circu mvented for nuclea r 
matter has still to be demonstrated. Miiller et al. (200C) 
have investigated whether nuclear matter can be formu- 
lated on a spatial lattice with nearest neighbor interac- 
tions, similar to the Hubbard model for high-T c super- 
conductors. Besides these theoretical efforts it is equally 
important to improve our understanding of the nuclear 
interaction in extremely neutron-rich nuclei (matter). 

Despite four decades of intense research, the astro- 
physical site of the r-process nucleosynthesis has yet not 
been identified. Recent astronomical and meteoric ev- 
idence points now to more than one source for the so- 
lar r-process nuclides, and it is clearly a major goal in 
the astrophysics community to solve this cosmic riddle. 
However, the puzzle will not be definitely solved if the 
nuclear uncertainties involved are not removed. This is 
even more necessary as recent research shows that the 
r-process is a dynamical process under changing astro- 
physical conditions. This implies dynamically changing 
r-process paths. To determine the paths one needs to 
know the masses of nuclei far-off stability accurately, be- 
sides the condition of the astrophysical environment. In 



a dynamical r-process, the astrophysical timescale will 
compete with the nuclear timescale, i.e. with the time 
needed for the matter flow from the seed nuclei to the 
heavier r-process nuclides. This nuclear timescale is set 
by the half-lives of the nuclei along the paths, in par- 
ticular by those of the longer-lived waiting point nuclei 
associated with the magic neutron numbers. The half- 
lives of waiting point nuclei are a very illustrative exam- 
ple for the need of reliable experimental data: Modern 
global theoretical models predict half-lives at the waiting 
points with a spread of nearly one order of magnitude 
and only data can decide. For the N = 50 and N = 82 
waiting point s such data exist for a few key nuclides (e.g. 



Pfciffer et al 



2001D , but not for the N = 126 nuclei. 
R-process nucleosynthesis as well as other astrophys- 
ical processes will tremendously benefit from future ex- 
perimental developments in nuclear physics. Worldwide 
radioactive ion-beam facilities, with important and dedi- 
cated programs in nuclear astrophysics, have just started 
operation or are under construction or in the proposal 
stage. These new facilities will boost our knowledge 
about nuclei far-off stability, they will determine astro- 
physically relevant nuclear input directly (e.g. masses and 
half-lives for the r-process, half-lives and cross sections 
for the rp- process, etc). But equally important, the ra- 
dioactive ion-beam facilities will guide and constrain the 
nuclear models, in this way indirectly contributing to 
the reduction of the nuclear inaccuracies in astrophysical 
models. Ultimatively nuclear physics can and will then 
become a stringent test and guidance for astrophysical 
theories and ideas. 
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